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ABSTRACT 

We examine the spheroid growth and star formation quenching experi¬ 
enced by galaxies since z ~ 3 by studying the evolution with redshift 
of the quiescent and spheroid-dominated fractions of galaxies from the 
CANDELS and GAMA surveys. We compare the observed fractions with 
predictions from a semi-analytic model which includes prescriptions for 
bulge growth and AGN feedback due to mergers and disk instabilities. We 
facilitate direct morphological comparison by converting our model bulge- 
to-total stellar mass ratios to Sersic indices. We then subdivide our popu¬ 
lation into the four quadrants of the sSFR-Sersic index plane and study the 
buildup of each of these subpopulations. We find that the fraction of star 
forming disks declines steadily, while the fraction of quiescent spheroids 
builds up over cosmic time. The fractions of star forming spheroids and 
quiescent disks are both non-negligible, and stay nearly constant over the 
period we have studied. Our model is qualitatively successful at reproduc¬ 
ing the evolution of the two “main” populations (star forming disks and 
quiescent spheroids), and approximately reproduces the relative fractions 
of all four types, but predicts a stronger decline in star forming spheroids, 
and increase in quiescent disks, than is seen in the observations. A model 
with an additional channel for bulge growth via disk instabilities agrees 
better overall with the observations than a model in which bulges can grow 
only through mergers. We also examine the relative importance of these 
different physical drivers of transformation (major and minor mergers and 
disk instabilities). 
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- galaxies: star formation 


1 INTRODUCTION 

The mechanisms by which galaxies are transformed 
and evolve over time, both in terms of their star for- 
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mat ion rates and their morphologies, are still not 
clearly known. At low redshift, the distribution of 
galaxy colors is bimodal ( Baldry et al.[ 2004; Bell et al. 


2004b). This division of galaxies into the star forming 
“blue cloud” and the quiescent “red sequence” can be 
observed most clearly in the color-magnitude and spe- 
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cific star formation rate (sSFR)-stellar mass planes 


(Baldry et al. 2004 Brinchmann et al. 2004 Kauff- 

mann et al. 2003 Strateva et al. 2001). Additionally, 
star forming galaxies can be said to occupy a “star 
forming main sequence,” a correlation between the 
star formation rate and the stellar mass of star form- 
ing galaxies (Noeske et al. 2007| Daddi et al. 2007s 

Elbaz et al. 2007| Rodighiero et al. 2011). Galaxies 

that are part of the red sequence have a wider range of 
star formation rates, although they do exhibit a corre¬ 
lation between mass (or luminosity) and color, where 
more massive galaxies tend to be redder (|Bernardi 

et al. 2003; Gallazzi et al. 2006; Peng et al. 2010; 

Brammer et al. 2011; Muzzin et al. |2013b). In ad- 

dition to the bimodality due to stellar populations, 
there is also a bimodality in the structure of galaxies 

(Kauffmann et al.|]2003). Often characterized by the 

bulge-to-total luminosity or mass ratio, or light pro¬ 
file parameterizations such as the Sersic index, galaxy 
morphology tends to be correlated with the star for¬ 
mation activity in the galaxy. Galaxy disks tend to be 
bluer than bulges (Peletier k Balcells|1996 Bell et al. 

2004a) and galaxies that are part of the blue cloud 
are more likely to be disk-dominated, while galaxies 
that are members of the red sequence are more likely 
to have more prominent bulges, or to have the con¬ 
centrated light profiles that are characteristic of early 
type galaxies (Blanton k Moustakas 2009; Schimi- 

novich et al. 2007 Bell|2008 Cheung et al.|2012). 

Large surveys have shed light on how the galaxy 
population evolved over a large fraction of the age of 
the universe. These observations have shown that the 
bimodality seen in the local universe is in place even 

at z ~ 2 — 3 (Brammer et al.|2009 2011 Muzzin et al. 

2013b). Analysis of the buildup of stellar populations 
from high redshift to the present reveals that the stel¬ 
lar mass contained in objects in the blue cloud has 
remained relatively constant, while the stellar mass 
represented by galaxies on the red sequence has grown 
significantly; this implies that blue star forming galax¬ 
ies are in fact being transformed into red, quiescent 

ones (Bell et al.|f2004b Borch et al. 2006 Bell et al. 

2007 |Faber et al. 2007). The mechanism responsi- 

ble for this “quenching” (or turning off of star for¬ 
mation) is not so clear. New information about the 
evolution of galaxy structure and morphologies has re¬ 
cently been gleaned from observations using the Hub¬ 
ble Space Telescope. Recent work suggests that quies¬ 
cence is intimately tied to the presence of a bulge com- 
ponent (Wuyts et al. 2011b Lang et al.| 2014| Bluck 

et al.|2014 McIntosh et al.|2014). Moreover, observa- 

tions have revealed a population of compact spheroid- 
dominated star forming galaxies at z ~ 2, which may 
be the progenitors of the quiescent, elliptical galax- 

ies we see today (Wuyts et al.|2011b; Whitaker et al. 

2012b; Barro et al. 2013; Williams et al.||2014| |Barro 


the mechanisms responsible for quenching, morpho¬ 
logical, and size evolution are connected. 

There have been several mechanisms proposed to 
explain galaxy quenching. One of the most popular 
scenarios involves feedback due to active galactic nu¬ 
clei (AGN). AGN feedback can be broadly divided 


into two regimes: the radiatively efficient “quasar” or 
“bright” mode, which is proposed to drive a powerful 
wind which expels gas from the galaxy, and the “ra¬ 
dio” or “maintenance” mode, which heats gas in the 
galactic halo, preventing it from cooling and form¬ 


ing stars (Somerville k Dave 2014 and references 


therein). This AGN activity can be driven either by 
galaxy mergers (Ellison et al.||2011 Silverman et al 
2011]) or in situ processes such as disk instabilities 
( Bournaud et al.|[2011 Dekel k Burkert||2014 ). Both 
of these processes lead to rapid transfer of angular mo¬ 
mentum and the growth of a bulge component. Virial 
shock heating is another proposed mechanism: during 
collapse, gas can be heated via the conversion of grav¬ 
itational potential energy into kinetic energy ( |White 
k Rees|[l978 ). Above a (redshift dependent) critical 


halo mass of ~ 1O 12 M 0 , this shock heating may be 
able to keep a substantial fraction of the halo gas hot, 
leading to quenching (Birnboim k Dekel| 2003 Keres 
etaL|2005t . 


While this does not seem directly related to the 
presence of a bulge component, it is clear from obser¬ 
vations that galaxies residing in halos above 1O 12 M 0 
are more likely to be bulge-dominated than disk- 
dominated ( Dekel et al.|2009 Woo et al.|2015 ). There 
is also the possibility that the presence of a (signifi¬ 
cant) bulge may itself stabilize the disk against local 
instabilities, thus making star formation less efficient, 


an effect known as morphological quenching (Mar 
tig et al. 2009). Finally, there is a suite of processes 


connected with dense environments, including tidal 
and ram pressure stripping and harrassment. These 
are often collectively referred to as “environmental 
quenching” ( Oemler|1974| Dressier 1 198 0; Balogh et al. 


2004 Tinker k Wetzel 2010 Peng et al. 2010), and 


they likely primarily affect satellites orbiting within 
a larger halo. These processes probably operate on a 
different timescale, and lead to different sorts of mor¬ 
phological transformation, than the ones described 
above. In this paper, we focus on field galaxy envi¬ 
ronments, so environmental processes are likely to be 


sub-dominant. See also McIntosh et al. (2014) for a 


summary of proposed quenching processes. 

If quenching and morphological transformation 
are (in most cases) intimately tied to each other, 
then galaxies which seem to be the “outliers” in this 
picture may be of particular interest: the quiescent 
disk-dominated galaxies and star forming spheroid- 
dominated galaxies. These populations are smaller 
than those of star forming disk-dominated and quies¬ 
cent spheroid-dominated galaxies, although they are 


not insignificant in size ( McGrath et al.|2008 van der 
Wei et al.|2011 ). Schawinski et al. ( 2014| ) did an anal¬ 


ysis of galaxies in the local universe that occupy the 
“green valley,” the region in between the blue cloud 
and the red sequence on the color-magnitude dia¬ 
gram, using observational data from SDSS (York et al. 


2000) and GALEX (Martin et al. 2005). They used 


morphology classifications from Galaxy Zoo (Lintott 
et al.|2008 2011) and determined that there were two 


distinct paths through the green valley, one taken 
by galaxies that leave the blue cloud as disk dom¬ 
inated systems, the other by galaxies that transi- 
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tion as bulge-dominated systems (see also the work 
of ( Smethurst et al.j|2015 ). The path associated with 
spheroid-dominated galaxies is consistent with work 
that suggests that bulge growth precedes quiescence 


(Wuyts et ah 2011b Wong et ah 


2012 


Lang et al. 


2014), while the path taken by disk-dominated galax¬ 


ies may explain the slowly growing population of qui¬ 
escent disk galaxies in the local universe which have 
been cut off from their gas supply but suffered noth¬ 
ing catastrophic to destroy or use up their existing 
gas reservoirs. Barro et al. (2013) and Woo et al. 


(2015) similarly identify different scenarios of quench¬ 


ing based on structural evolution. McIn tosh et al.| 
(2014) also describe a path associated with spheroid- 


dominated galaxies when examining what they call 
“recently quenched ellipticals” in the Sloan Digital 
Sky Survey. In this light, these “outlying” populations 
are possibly much more important to the overall pic¬ 
ture than originally suspected, and may be an indi¬ 
cation that quenching is caused by multiple physical 
processes. 

Aiding in the investigation of galaxy formation 
and evolution are state of the art numerical simu¬ 
lations and semi-analytic models. N-body dark mat¬ 


ter simulations such as Bolshoi (Klypin et al. 2011) 


are an invaluable tool when testing the predictions 
of our currently favored Lambda Cold Dark Matter 
cosmological model, ACDM. Semi-analytic models, or 
SAMs, plant galaxies in merger trees assembled in 
dark matter N-body simulations or constructed us¬ 
ing techniques based on the Extended Press Schechter 
formalism. By following the evolution of these galax¬ 
ies within the backbone of the dark matter history, 
accounting for physical processes such as gas accre¬ 
tion and cooling, star formation, merging and feed¬ 
back with physically motivated recipes, population 
statistics for a cosmological sample of galaxies can be 
generated quickly and with minimal computational 


resources (jKauffmann et al. 

19931 |Cole et al.l 

1994 

Somerville & Primack 

1999 

De Lucia et al.! 

2006 


et al. 


2012 


Somerville et alj 2008a[ |Guo et ah |2011| |Somerville| 


Porter et al. 2014a). SAMs have been 


used to study the evolution of star formation and the 
buildup of spheroid-dominated galaxies and to gauge 
which processes are especially impor tant. |De Lucia| 
et al. (2006) and Benson & Devereux ( 20101) have in¬ 


vestigated the buildup of spheroid-dominated galaxies 
in SAMs with cosmic time. These analyses and oth¬ 
ers have reiterated that bulge growth often appears 
to be connected to the cessation of star formation 
and also demonstrate that the two main channels for 
bulge growth are mergers and disk instabilities, both 
of which appear to be important, although their de¬ 
gree of importance may change with redshift, galaxy 
mass and environment, and is also model dependent 
(Parry et al. 2009 De Lucia et al.|20~ Fontanot et al. 


2012 Porter et ai.||2014a ). 


Recently, Porter et al. (2014a) compared the pre¬ 
dictions of the latest version of the “Santa Cruz” SAM 


(Somerville & Primack 1999; Somerville et al.f2 008a) 
for the z = 0 stellar mass function divided by mor¬ 
phology with available observations, and found fairly 
good agreement. They found that adding a prescrip¬ 


tion for bulge growth via disk instabilities brought the 
model into better agreement with the observed galaxy 
stellar mass function of spheroid-dominated galaxies 
at intermediate masses (10.5 < log(M*/M 0 ) < 11.5). 
In addition, Porter et al. (2014a) developed a new 


model for predicting the radial sizes and velocity dis¬ 
persions of bulges formed via mergers or disk instabili¬ 
ties, based on a simple analytic model calibrated using 
numerical hydrodynamic simulations of binary galaxy 
mergers. Their model reproduces the observed size- 
mass relation for spheroids and disks, and the evolu¬ 
tion of this relation from z ~ 2 to the present day (see 
also Somerville et al. in prep). |Porter et al.| (2014b) 
investigated the predictions of the same models for 
the correlation of the age and metallicity of stars in 
local spheroid-dominated galaxies with structural pa¬ 
rameters such as size and velocity dispersion. They 
found a strong correlation between both stellar pop¬ 
ulation parameters (age and metallicity) and internal 
velocity dispersion, in agreement with observations. 
They found no correlation between age and radius, 
and a weak correlation between metallicity and ra¬ 
dius, also in agreement with observations of nearby 
early type galaxies. In this paper we follow up on the 
work by Porter and collaborators by directly study¬ 
ing the build-up of the spheroid-dominated popula¬ 
tion over cosmic time, and comparing with observa¬ 
tions of high-redshift galaxies. 

The conclusions of previous studies in the lit¬ 
erature regarding the spheroid-dominated fraction 
of galaxies in SAMs and the agreement with ob¬ 
servations are difficult to synthesize, because differ¬ 
ent analyses use different criteria to define spheroid- 
dominated galaxies both in the models and in the ob¬ 
servations. The bulge-to-total mass ratio (. B/T ) that 
is readily predicted in SAMs is difficult to measure ob- 
servationally. Therefore it has been difficult to make a 
direct comparison between model predictions and ob¬ 
servations previously. One of the important new fea¬ 
tures of this study is that we extend our models to 
predict a morphological quantifier that can be com¬ 
pared more directly with observations. We describe 
our new method in detail below. 

In this paper, we present new results quantifying 
the evolution of quenching and spheroid growth in ob¬ 
servations from z ~ 3 to the present, and also present 
new predictions of these same quantities from state- 
of-the-art semi-analytic models. We split galaxies ac¬ 
cording to their star formation rates and morpholo¬ 
gies and examine the buildup of the quiescent and 
spheroid-dominated fractions of galaxies. We then go 
further than studies in the past by subdividing into 
four populations: star forming disk-dominated galax¬ 
ies (SFDs), star forming spheroid-dominated galax¬ 
ies (SFSs), quiescent disk-dominated galaxies (QDs) 
and quiescent spheroid-dominated galaxies (QSs). We 
examine the evolution of the fraction of galaxies in 
each of these populations. Our low redshift observa¬ 
tional data (z rsj 0.06) come from the Galaxy and 


Mass Assembly Survey (GAMA Driver et al. 2009) 


and our higher redshift data (0.5 < z < 3.0) come 
from the Cosmic Assembly Near-infrared Deep Ex- 
tragalactic Legacy Survey (CANDELS; Grogin et al. 
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2011 Koekemoer et al.|2011 ). We use the Santa Cruz 
SAM of Somerville et al. (2008a) with updates as de¬ 


scribed in Somerville et al. (2012) and Porter et al 
(2014aJ. Our semi-analytic model includes the effects 
of AGN feedback and bulge growth triggered by merg¬ 
ers and (optional) disk instabilities. Another way in 
which our study is unique is the way we characterize 
the morphologies of our model galaxies: we convert 
our model output, bulge-to-total mass ratio, to Ser- 
sic index as described in Section 2 and Appendix A 
in order to facilitate a more direct comparison be¬ 
tween model and observed galaxy morphologies than 
has been carried out before. We will also examine in 
detail the histories of galaxies selected from each pop¬ 
ulation in order to shed light on the individual tracks 
that different types of galaxies move along as they 
evolve. The structure of the paper is as follows. In Sec¬ 
tion 2 we give an overview of our semi-analytic model 
and of the data sets with which we are comparing. In 
Section 3 we present a comparison of the evolution of 
these populations in the model and the observations. 
We present our discussion, in part informed by study¬ 
ing individual evolutionary tracks of galaxies from the 
model, in Section 4 and our summary and conclusions 
in Section 5. 


2 SEMI-ANALYTIC MODEL AND 

OBSERVATIONAL DATA 

2.1 The Semi-Analytic Model 

The SAMs used in this paper were first presented in 
Somerville & Primack (1999) and Somerville et al. 
( 2001)) , and significantly updated in |Somerville et ah 
( 2008a[ S08), |Somerville et~ah ( 2012[ S12) and Porter 
et al. ( 2014a[ PI4). The model includes prescriptions 
for the following physical processes: the hierarchical 
growth of structure in the form of dark matter merger 
trees, the heating and cooling of gas, star formation 
as governed by the empirical Kennicut-Schmidt law, 
the evolution of stellar populations, supernova feed¬ 
back, chemical evolution of the ISM and ICM due to 
supernovae, AGN feedback, and starbursts and mor¬ 
phological transformation due to galaxy mergers and 
disk instabilities. Here we will briefly summarize these 
processes, focusing mainly on the processes relevant 
to the evolution of star formation and morphology. 
For a more in depth description of the model, see 
S08 and P14. We assume a standard ACDM cos¬ 
mology (9 m — 0.27, Qa — 0.73, h — 0.7 ) and a 
Cha brier] (2003) initial mass function. Our adopted 
baryon fraction is 0.1658. Our cosmology was chosen 
to match that adopted by the Bolshoi simulation (de¬ 
tailed below) and is consistent with the Wilkinson Mi¬ 
crowave Anisotropy Probe (WMAP) 5/7-year results 
(Komatsu et al. l2059l[20lT] ). 

In this work, we use the CANDELS lightcones 
(Somerville et al. in prep) extracted from the Bolshoi 
dark-matter only N-body simulation (Klypin et al. 
2011 Trujillo-Gomez et al. |2011 ). Dark matter ha¬ 
los are identified using the ROCKSTAR algorithm 


has a force resolution and mass resolution of l/i -1 
kpc and 1.9 x 10 8 M Q , respectively. Merger trees are 
constructed for each halo in the lightcone using the 


method of Somerville & Kolatt (1999). There is no 


appreciable difference in results when using merger 
trees extracted from the N-body simulation (as done 
in Porter et al. 2014) as opposed to EPS (as we do 
here). For our lowest redshift bin, the lightcones rep¬ 
resent a very small volume so we simply use a low-z 
snapshot from the Bolshoi volume. 

When dark matter haloes merge, the central 
galaxy of the largest progenitor becomes the new cen¬ 
tral galaxy, while all other galaxies become satellites. 
Satellite galaxies are able to spiral in and merge with 
the central galaxy, losing angular momentum to dy¬ 
namical friction as they orbit. The merger time-scale 
is estimated using a variant of the Chandrasekhar for¬ 


mula from Boylan-Kolchin et al. (2008). Tidal strip¬ 


ping and destruction of satellites as described in S08 
are also included. 

Before the universe is reionized, each halo has a 
hot gas mass equal to the virial mass of the halo times 
the universal baryon fraction. The collapse of gas into 
low-mass haloes is suppressed after reionization due 
to the photoionizing background. We assume the uni¬ 
verse is fully reionized by z = 11 and use the results 


of Gnedin (2000) and Kravtsov et al. (2004) to model 


the fraction of baryons that can collapse into haloes of 
a given mass following reionization. Due to the galaxy 
mass range selected in this work, we do not expect our 
results to be sensitive to this prescription. 

When dark matter haloes collapse or are involved 
in a merger that at least doubles the mass of the 
progenitors, the hot gas is shock-heated to the virial 
temperature of the new halo. The rate at which this 
gas can cool is determined by a simple spherical cool¬ 
ing flow model. Assuming a monotonically decreas¬ 
ing density profile for the gas, and that denser gas 
cools faster, we can define a “cooling radius”, within 
which all gas is able to cool within some time t c 00 i, 
which we have defined as the halo dynamical time. 
The initial density profile is assumed to be that of a 
singular isothermal sphere, and the cooling radius is 
found by using the atomic cooling curves of |Suther^| 
land &; Dopita (1993). The cooling radius may be 
larger or smaller than the virial radius of the halo; 
when the cooling radius is larger, the cooling rate is 
limited only by the rate at which gas is infalling. The 
transition from r coo i > rvn to r co oi < rVir is associ¬ 
ated with the transition from “cold flows”, where cold 
gas streams into the halo along dense filaments with¬ 
out being heated, to “hot flows”, where gas is shock 
heated on its way in, forming a diffuse hot gas halo 
before cooling (Birnboim &; Dekel|2003 Dekel & Birn- 


of Behroozi et al. (2013a). The Bolshoi simulation is 
complete down to halos with Wire = 50 km/s, and 


boim||2006 Keres et al.|2005j . Note that in this way, 

virial shock heating (sometimes referred to as ‘halo 
mass quenching’) is included in our SAMs. However, 
it has been shown by many studies (both numerical 
and semi-analytic) that this effect alone is insufficient 
to create the observed population of massive quies¬ 
cent galaxies ( |Somerville fe D ave|20l4 , and references 
therein). 

Newly cooled gas collapses to form a rotationally 
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supported disk, the scale radius of which is estimated 
based on the initial angular momentum of the gas and 
the profile of the halo. We assume that angular mo¬ 
mentum is conserved and that the self-gravity of the 
collapsing baryons causes the inner part of the halo to 


contract ( 

Blumenthal et al. 1986 

Flores et al.||1993| 

Mo et al. 

1998). This method was shown to repro- 

duce the observed size-stellar mass relations of disks 

out to z r 

^ 2 in Somerville et al. (2008b). Spheroids 


can be created by mergers or disk instabilities. The 
sizes of spheroids formed in mergers are determined 
by the stellar masses, sizes and gas fractions of the two 
progenitors, as described in P14. The size of spheroids 
formed in disk instabilities is determined by assuming 
that they form from the center of the exponential stel¬ 
lar disk; the radius is simply the radius that contains 
the amount of mass that is to be transferred from the 
disk to the bulge (again, see P14 for details). 

There are two modes of star formation in the 
model: a “normal” mode that occurs in isolated disks 
and a “starburst” mode that occurs as a result of a 
merger or internal disk instability, which will be dis¬ 
cussed in more depth below. The normal mode follows 
the Schmidt-Kennicutt relation (Kenn icutt|1998 ) and 
assumes that gas must be above some fixed critical 
surface density (the adopted value here is 6M 0 /pc 2 ) 
in order to form stars. 

Exploding supernovae and massive stars are ca¬ 
pable of depositing energy into the ISM, which can 
drive outflows of cold gas from the galaxy. We assume 
that the mass outflow rate is proportional to the SFR 
and decreases with increasing galaxy circular veloc¬ 
ity, in accordance with the theory of “energy-driven” 
winds. Some ejected gas is removed from the halo 
completely, while some is deposited into the hot gas 
reservoir of the halo and is eligible to cool again. The 
gas that is driven from the halo entirely is combined 
with the gas that has been prevented from cooling 
by the photoionizing background and may later reac¬ 
crete back into the halo. The fraction of gas which is 
retained by the halo versus the amount that is ejected 
is a function of halo circular velocity as decribed in 
S08. 


Heavy elements are produced by each generation 
of stars, and chemical enrichment is modelled simply 
using the instantaneous recycling approximation. For 
each parcel of new stars dm*, a mass of metals dMz = 
ydm * is also created, which is immediately mixed with 
the cold gas in the disk. The yield y is assumed to be 
constant and is treated as a free parameter. Supernova 
driven winds act to remove some of this enriched gas, 
depositing a portion of the created metals into the hot 
gas or outside of the halo. 


2.1.1 Mergers and Starbursts 

Mergers between galaxies are assumed to remove an¬ 
gular momentum from stars and gas in the disk 
and drive material towards the center, building up a 
spheroidal component. In our model, this spheroidal 
component is formed instantaneously. In principle 
this could affect our results by forming bulges more 
quickly than they should form in the real universe. 


However, actual bulge formation time scales (~tdyn) 
are quite short compared with the times associated 
with our redshift bins, so we don’t expect this to have 
much of an effect. 

Mergers also trigger a starburst, the efficiency 
of which depends on the gas fraction of the central 
galaxy and the mass ratio of the two progenitors. The 
time scale of the burst is also determined by properties 
of the progenitor galaxies. The parameterization is 
based on hydro dynamical simulations of binary merg¬ 
ers between disks (Hop kins et al.|2009b ). Simulations 
show that the closer the mass ratio of the progenitors 
is to one (or how “major” the merger is) and the more 
gas-poor the merger is, the more efficient it is at re¬ 
moving angular momentum from the gas and driving 
it into the nucleus, and scattering disk stars into a hot 
spheroid component ( Cox et al.|2006 Robertson et al. 


2006). The gas fraction dependence can be understood 


as follows: if the progenitors are very gas-rich, there 
is not enough stellar mass to create a torque on the 
gas, making it difficult for the gas to shed angular mo¬ 
mentum and collapse inward (Hopkins et al. 2009b). 
S08 parameterized the burst efficiency only as a func¬ 
tion of mass ratio, but S12 and P14 introduced the 
gas fraction dependence in accordance with Hopkins 


et al. (2009a). Stars that are formed as part of the 


starburst are added to the spheroidal component, as 
are 80% of the stars from the merging satellite galaxy. 
The other 20% are deposited into a diffuse stellar halo 
component. 


2.1.2 Disk Instabilities 


Disk material can also be converted into a spheroidal 
component as a result of internal gravitational insta¬ 
bilities. A pure disk without a dark matter halo is very 
unstable to the formation of a bar or bulge, while mas¬ 
sive dark matter haloes tend to stabilize a thin, cold 
galactic disk (Ostriker & Peebles 1973| Fall &; Efs- 
tathiou|1980 ). When the ratio of dark matter mass to 
disk mass falls below a critical value, the disk can no 
longer support itself and material collapses into the 


inner regions of the galaxy (Efstathiou et al. 1982). 


Here we adopt an avenue for bulge growth due to disk 
instability, based on a Toomre-like stability criterion. 


Following Efstathiou et al. (1982), Mo et al. (1998), 
P14 and many other works, we define the stability 
parameter as 




_ * max ( 1 ^ 

£disk “ (GM disk /r disk )V3 W 

where E ma x is the maximum circular velocity of 
the halo (used as a proxy for the maximum circular 
velocity of the disk), raisk is the scale length of the 
stellar disk and Mdisk is the stellar mass of the disk. 
This is identical to the “Stars DI” disk instability cri¬ 
terion introduced in P14. Whenever eaisk < e cr it, the 
disk is considered to be unstable. The value of e cr it in 
numerical simulations of isolated disks has been found 
to be in the range of 0.6 — 1.1, with disks containing 
stars and cold gas having a lower threshold than pure 
stellar disks ( Efstathiou et al.||l982 Mo et al!||1998 ). 
We set e cr it = 0.75 as in P14, where this value was 
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chosen to match the observed fraction of spheroid- 
dominated galaxies at z — 0. When the disk becomes 
unstable, stellar mass is moved from the disk to the 
bulge until eaisk = e cr it- The gas in the disk is not af¬ 
fected. The “Stars+Gas DI” model of P14 included 
gas in determining the stability of the disk and also 
moved some gas to the bulge component to feed the 
central supermassive black hole when the disk became 
unstable. However, the results for the two approaches 
were very similar. Again, the creation of the bulge 
component is instantaneous. While we are aware that 
this implementation of disk instability is crude and 
perhaps does not capture all of the relevant physics, 
this is an approach that is commonly used in the lit¬ 
erature. One of the goals of this work is to explore 
how important bulge growth through disk instabili¬ 
ties might be, in order to guide future investigations. 
Later, we discuss more physical models of disk insta¬ 
bility and how including them might affect the results 
of this study. 

It is worth noting that we also do not account for 
the possibility that a previously existing bulge may 
help stabilize the disk against another instability. Be¬ 
cause of this (and the fact that we only move as much 
material as needed to restabilize the bulge) it is pos¬ 
sible (even common) for disks to develop chronic in¬ 
stabilities which lead to the steady growth of a bulge 
component. 

Below we present our results for versions of 
the SAM both with the disk instability prescription 
turned on (DI model) and off (noDI model). The DI 
model is our fiducial model, however, and unless oth¬ 
erwise noted, it is the DI model that is shown. 


2.1.3 Black Hole Accretion and Feedback 


Galaxies ar e initially seeded with a massive black hole 
of 10 4 Mq ( Hirschmann et al.|2012 ). When two galax¬ 
ies merge as described above, their central black holes 
are assumed to merge as well, after which the new 
central black hole of the merger remnant engages in a 
bout of feeding and radiatively efficient, or “quasar” 
mode, AGN activity. During this time, the black hole 
accretes at its Eddington limit. As the black hole ac¬ 
cretes and radiates, it deposits energy into the sur¬ 
rounding medium until it reaches a critical mass which 
corresponds to the energy which would stop accre¬ 
tion and begin driving an outflow, such as those seen 


in many recently merged systems) Rupke & Veilleux 
Emo nts et al.pOM ). The black hole effectively 


2013 


starves itself of material, as its accretion rate declines 
as a power law, in accordance with the results of |Hop-| 
kins et al. (2006). We follow the hydrodynamical bi¬ 
nary merger simulations of Hopkins et al. ( 2007| for 
our definition of the critical mass, M cr it, at which the 
black hole accretion rate enters the declining phase, 
and Affinal ? at which the black hole stops feeding. If 
the newly merged black hole is already more massive 
than Mfi na i, there is no accretion event. We note that 
our predicted final black hole and bulge masses are 
consistent with the observed Mbh — Afbuige relation 


( Somerville et al.|2008a Hirschmann et al. 2012). 

A bout of black hole accretion and AGN activ¬ 


ity can also be triggered by a disk instability. When 
disk mass is transferred to the bulge as previously 
described, we assume the black hole accretes a gas 
mass equivalent to some fraction of that mass. Fol¬ 
lowing Hirschmann et al. (2012), we set this term to 
be /fuei,Di = 0.002, which leads to good agreement 
with the observed number density of low-luminosity 
AGN. The black hole can continue to accrete until 
this fuel is consumed. 

The black hole is also able to feed and effect feed¬ 
back in the “radio” or “maintenance” mode. In this 
mode the black hole feeds via Bondi-Hoyle accretion 
from the hot halo ( Bondi|1 952). The accretion is usu¬ 
ally significantly sub-Eddington. This feedback mode 
is associated with giant radio jets which heat the sur¬ 
rounding gas, preventing it from cooling and forming 
stars. Once the accretion rate is determined, a cou¬ 
pling constant determines how effectively the energy 
released couples to the surrounding gas. The radio 
mode heating rate is then calculated and subtracted 
from the cooling rate described above. 


2.1.4 Computing Sersic Indices 


This work involves comparing the morphologies of 
model galaxies with Af* > 10 lo Af Q with those of ob¬ 
served galaxies. From the SAM, we can easily cal¬ 
culate the bulge-to-total stellar mass ratio, and we 
have predictions for the radii of the galactic bulge 
and disk. There are many different methods used to 
classify the morphologies of observed galaxies. One 
commonly used method is to fit the light profile with 
a Sersic function, resulting in the determination of the 
Sersic index ( Sersic|19 63). In an effort to put the ob¬ 
servations and model on equal footing, we have con¬ 
verted our model outputs to a Sersic index using a 
lookup table which takes in the bulge-to-total mass 
ratio and bulge radius to disk radius ratio and gives 
an effective radius and Sersic index. This lookup table 
was generated by fitting Sersic indices and effective 
radii to synthetic bulge+disk systems (n— 1 for disks 
and n— 4 for bulges) for a range of different bulge-to- 
total mass ratios and bulge radius to disk radius size 
ratios. The values that come out of the lookup table 
are discrete for obvious reasons, so we use a 2D in¬ 
terpolation of the table to generate our Sersic indices 
and effective radii. More information and some tests 
of our approach can be found in Appendix A. 


2.2 Observational Samples 

2.2.1 CANDELS 


In this work, we make use of Hubble Space Tele¬ 
scope/Wide Field Camera 3 (HST/WFC3) observa¬ 
tions of galaxies taken as part of the Cosmic Assem¬ 
bly Near-infrared Deep Extragalactic Legacy Survey 
(CANDELS; |Grogin et al. 2011 |Koekemoer et aT 


[ 2 OII). These observations span two of the five CAN- 


DELS fields: the Ultra-Deep Survey (UDS; [Lawrence 


et al.|2 007) and the Great Observatories Origins Deep 


Survey South (GOODS-S; Giavalisco et al. 2004). 
With the combined strengths of galaxy selection in the 
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F160W (H) band and the availability of rich multi¬ 
wavelength datasets, the CANDELS catalogs afford 
us an unprecedented study of galactic structure and 
the rise and fall of star formation activity from high 
redshift toward the present day. 

Here, we give only an outline of the galaxy cata¬ 
logs used for our analysiiQ The galaxies were drawn 
from the CANDELS catalogs for UDS and GOODS-S; 
for more details we refer the reader to lGalametz et all 


dust extinction assuming the Calzetti law (Calzetti 


(2013) and Guo et al. (2013), respectively. Briefly. 


the UV-to-NIR multiwavelength photometric cata¬ 
logs were computed using the template-fitting method 
TFIT ( Lee et ah]|2012 Laidler et al.||2007 ) which al¬ 
lows us to consistently merge datasets with signifi¬ 
cantly different spatial resolution. 

Photometric redshifts were determined following 


the method described in Dahlen et al. (2013) which 


combines redshift probability distributions from sev¬ 
eral different codes using a Bayesian approach to im¬ 
prove the precision and reduce the number of catas¬ 
trophic outliers. For a sample of 480 spectroscopi¬ 
cally confirmed GOODS-S galaxies without an active 
galactic nucleus, the median and standard deviation 
of (Zphot - 2 sp ec)/(l + ^spec) are approximately -0.014 
and 0.045, respectively. 

The stellar masses were drawn from the catalog 


presented in Santini et al. (2015). The catalog includes 


stellar masses computed using different SED fitting 
codes and modeling assumptions. For this work, we 


adopted stellar masses computed using FAST (Kriek 


et al. 2009) with Bruzual & Chariot (2003) stel 


lar population synthesis models, a Chabrier initial 


mass function Chabrier (2003), exponentially declin¬ 


ing star formation histories, solar metallicity and the 
Calzetti dust extinction law ( Calzetti] 2001). The total 
star formation rate for galaxies detected in the mid- 
and/or far-infrared is defined as SFRtotai=SFRuv+m 
= SFRuv+SFRir, where SFRuv = SFR2800 is the 
unobscured (and therefore uncorrected for dust ex¬ 
tinction) component derived from L 2800 (the lumi¬ 
nosity at 2800A), and SFRir is the integrated Lir- 
based obscured component. The integrated (or total) 
infrared luminosity, Lir, is itself derived from fitting 
24 /im fluxes using Chary & Elbaz (2001) templates 
and a calibration determined from Herschel data (see 
El baz et al.|2011 ). Two major assumptions underlying 
this mapping from 24 /am flux to the total IR (8/xm to 
1000/xm) luminosity are (1) the IR SEDs of galaxies 
do not evolve significantly with redshift, and (2) emis¬ 
sion from dust heated by an obscured AGN does not 
significantly increase the 24 /xm fluxes (again, consult 


tions). 

To be more 


Elbaz et al. (2011) for the validity of these assump- 


10 (Lir+3.3L280o), based on 


concrete, SF Ruv + ir = 1.0 9 x 


Kennicutt 


(1998) and 


B ell et al.| p 005|, and assuming a |Chabrier| (2003) 


IMF (see Barro et al.| ( |2011 ) for more information). 
For galaxies without an infrared detection and thus no 
total Lir estimate, we instead corrected SFRuv for 


1 All CANDELS catalogs are avail¬ 
able at the Rainbow Database- 

http:/ /arcoiris.ucolick.org/Rainbow_navigator_public/ 


2001), giving us a comparable estimate of SFR to tai- 
SFRuv,corr = SFRuv • 10°' 4 ' * 1 - 8 ' Av . The optical extinc¬ 
tion parameter Ay for each galaxy was determined 
by FAST, and the factor of 1.8 in the exponent is 
the attenuation parameter, k(A = 2800). We refer the 
reader to Wuyts et al. ( 2011a| and |Perez-Go nzalez 
et al. (2008) for further details about the derivation 


of total star formation rates, including a wealth of 
comparisons between different estimates of SFR to tai- 
The structural measurements of the observed 
galaxies (Sersic indices, in particular), as computed 
by GALFIT ( |Peng et al.| 2002] |van der Wei et al. 
2012} from the HST/WFC3 F160W (H-band) images, 
are used throughout this work. To secure the robust¬ 
ness and completeness of our sample, we consider only 
those galaxies with ?tifi 60 W < 25, M* > 10 lo M©, and 
GALFIT flag equal to 0 (good fits only). For the red¬ 
shift range considered in this paper (0.5 < z < 3.0), 
these three major selection cuts leave us with the 
following numbers of galaxies in each field: 1123 in 
GOODS-S and 1594 in UDS. The GALFIT high- 
quality flag cut was applied last; our sample sizes be¬ 
fore that particular cut were 1333 and 1798 galaxies 


in GOODS-S and UDS, respectively. Guo et al. (2013) 
and Galametz et al. (2013) respectively show that the 


GOODS-S and UDS samples are complete after tak¬ 


ing into account our mFi 60 \v cut. Guo et al. ([2013) 


further show that the completeness in GOODS-S also 
depends on morphology by splitting a synthetic com¬ 
parison sample into disk-dominated and spheroid- 
dominated subsamples (based on the Sersic index). 
Nevertheless, our morphology-dependent subsamples 
should continue to remain complete given our mFmow 
and mass selection criteria. 


2.2.2 GAMA 

Given the small volume probed by CANDELS at 
z < 0.5, and the need to compare our high-redshift 
results to those obtained from a robust low-redshift 
anchor point, we incorporate multi-wavelength data 
from the Galaxy And Mass Assembly (GAMA) sur¬ 
vey into our analysis. GAMA is a large (144 deg 2 ) 
spectroscopic survey that builds on the legacy of the 
Sloan Digital Sky Survey (SPSS;[Yor k et aL[2C )00) and 
the Two-degree Field Galaxy Redshift Survey (2dF- 
GRS; Colless et al. 2001| , reaching a limiting mag¬ 
nitude of r < 19.8 mag with > 98% spectroscopic 
completeness (see Driver et al. ( |2011 ) for a review of 
the first three years of GAMA). It is this intermedi¬ 
ate depth and high spectroscopic completeness, com¬ 
bined with photometry spanning a large wavelength 


range (1 nm to lm Liske et al. 2014)), in contrast 
to wider and shallower spectroscopic surveys in the 
past, that makes GAMA a unique survey, and an ideal 
complement to the CANDELS dataset. In addition 
to conducting its own spectroscopic observations with 
the AAOmega spectrograph on the Anglo-Australian 
Telescope, GAMA has assembled existing and is pur¬ 
suing new spectroscopic and imaging data in collab¬ 
oration with several other independent surveys (see 
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Baldry et al. (2010), and Tables 1 and 4 of Liske et al. 


(2014)). 

Specifically, this work makes use of Data Re¬ 
lease 2 (DR2) from the GAMA survey. DR2 pro¬ 
vides, among other things, local bulk flow-corrected 
redshifts (see Section 2.3 of Baldry et al.||2012 ), stel¬ 
lar masses (Taylor et al. 2011), Ha-based star for¬ 


mation rates (Gunawardhana et al. 2013 Hopkins 
et al.|2 013), observed and rest-frame photometry (Hill 


et al.||2011), and GALFIT structural measurements 
we use r-band fits). These data 


(Kelvin et al. 


2012 


are released for 72,225 objects distributed over three 
GAMA regions: two 48 deg 2 fields with limiting mag¬ 
nitude r < 19.0 mag, and one 48 deg 2 field with lim¬ 
iting magnitude r < 19 A ma g, giving a total survey 


volume of 144 deg 2 . 


Liske et al. 


(2014) provide an 


overview of DR2 for GAMA as well as further infor¬ 
mation about the survey’s progress. 


For the sake of consistency, the Ha-based SFR 
for each GAMA galaxy, SFRhcx, was converted from 


a Salpeter IMF basis to a Chabrier IMF basis(Bell 


et al. 2005). The SFRna measurements provided by 


GAMA should not, in principle, deviate greatly from 
the SFRuv+ir and SFRuv,corr prescriptions that we 
have adopted for CANDELS galaxies because the Ha 
luminosity (Ln a ) estimates for GAMA galaxies have 
been corrected for dust extinction and Balmer stel¬ 
lar absorption (Guna wardhana et al.|[2Q13 ). Further¬ 
more, we applied aperture corrections to the stellar 
mass and rest-frame magnitude estimates to account 
for the fraction of mass (or flux) that falls outside 
of the r-band (SExtractor AUTO |Ber tin & Arnouts 


1996) aperture used for aperture-matched photome¬ 
try (Liske et al. Pont . 


In order to ensure robustness and consistency, we 
apply the following selection cuts to the GAMA DR2 
catalog: M* > 10 10 M©, r-band GALFIT flag equal 
to 0 (good fits only), and high-quality redshifts only 
(“NQ”> 3 in the DR2 catalog). As mentioned ear¬ 
lier, the limiting magnitudes are r = 19.0 mag for 
two fields, and r = 19.4 mag for the third field. We 
consider only galaxies with local bulk flow-corrected 
redshifts between 0.005 and 0.12; the lower limit is to 
prevent stellar contamination in our galaxy sample (in 
fact, the lowest redshift in our sample after applying 
the aforementioned cuts is z ~ 0.0059), and the latter 
limit is the maximum redshift at which we are com¬ 
plete for M* ~ 10 lo M© galaxies. These selection cuts 
leave us with a total of 5112 GAMA DR2 galaxies. 
Before the GALFIT high-quality flag cut, but after 
all other cuts, our GAMA sample size was 5977. 


For each galaxy in our final GAMA sample, 
we derive completeness correction weights using the 


Usurvey/Umax weighting technique ( 

Schmidt|l968 and 

see Section 3 of Taylor et al. 2015 

) • ^survey is the to- 


tal comoving volume contained within the 144 square 
degree GAMA survey (taking into account the sur¬ 
vey geometry, and considering only our chosen red¬ 
shift slice, 0.005 < z < 0.12). U ma x is the maximum 
comoving volume within which a given galaxy could 
have been detected, again taking into account the sur¬ 
vey geometry (i.e., using the provided z ma x, 19.0 and 
£max,i 9.4 for each galaxy, corresponding to the three 


48 square degree GAMA fields with different r-band 
selection limits). Uurvey/Umax equals 1 for galaxies 
for which we are complete, and it is greater than 1 
for galaxies for which we are incomplete. The me¬ 
dian weight is 1.0, and only 187 galaxies (~3%) have 
a weight greater than 1.0 (with the maximum value 
being about 35). As expected, our GAMA sample is 
~97% complete. 


3 RESULTS 

We now examine how well our model (with the disk in¬ 
stability prescription turned both on and off) matches 
the buildup of the quiescent and spheroid-dominated 
fraction of observed galaxies. We then subdivide the 
model and observed populations further and exam¬ 
ine the buildup of the four quadrants of the sSFR- 
Sersic index plane. In this way we can assess where 
our model is succeeding and failing in transforming 
galaxies in terms of their star formation rates and 
morphologies. 


3.1 Quiescent Fraction 

3.1.1 Dividing by sSFR 


Our first step is to split galaxies into star forming 
and quiescent populations. We preferred not to sim¬ 
ply divide our population by eye and sought an au¬ 
tomated process which would divide our galaxies in 
each redshift bin in a reasonable way. One approach 
used in the literature is to divide at a specific star 
formation rate sSFR(z) = l/[3tn(z)], where tn(z) is 
the Hubble time at the redshift of interest. This di¬ 
vider in sSFR is roughly equivalent to the division of 
galaxies into star forming and quiescent on the UVJ 
color-color diagram as described in Whitaker et al.| 
(2012a) and Muzzin et al. (2013b). This division line, 
as well as others that we attempted to use, all shared 
the same problem: the distribution of sSFRs in the 
model and from observations is somewhat different, 
especially at z>2.2, so dividing lines which made a 
reasonable cut for model galaxies did not work as 
well for observed galaxies and vice versa. The sSFR 
distribution of model galaxies is not as bimodal as 
it is for the observations; rather than having a sec¬ 
ond peak at very low sSFR, our model distribution 
tails off. We don’t expect this to significantly affect 
the results of this work as the star formation rates in 
question are already very low (our model galaxies are 
being quenched; their sSFRs just aren’t distributed in 
quite the same way as the observations) and any new 
stars formed shouldn’t change the structural parame¬ 
ters with which we concern ourselves later. Still, this 
makes defining quiescence by examining the trough 
between populations somewhat difficult. We instead 
seek to define our dividing line in relation to the star 
forming main sequence, which leads us to a slightly 
different issue. 

The star formation rates of the observed galax¬ 
ies are systematically slightly higher than those of 
the model galaxies, so a typical observed star forming 
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galaxy (one which we would say occupies the main se¬ 
quence of star formation) has a different sSFR than a 
corresponding model galaxy. The dependence of sSFR 
on stellar mass for star forming galaxies is also steeper 
for observed galaxies than for model galaxies. This 
may point to a deficiency in some of our prescriptions 
for star formation and/or stellar feedback (see the 
discussion in |Somerville &; Dave||2014[ and references 
therein). However, in this work we are concerned with 
broadly distinguishing between star forming and qui¬ 
escent, and with the processes responsible for moving 
galaxies fairly dramatically off of the main sequence. 
Therefore as long as we define our dividing line rel¬ 
ative to the main sequence in the models and in the 
observations, our analysis should be robust. 

To deal with these issues, we introduce a method 
to calculate a dividing line between star forming and 
quiescent galaxies which we apply to both the obser¬ 
vations and the model galaxies; however, the actual 
normalization, slope, and redshift dependence of the 
dividing line are not the same for the model and the 
observational samples. Geha et al. (2012) has shown 
that, in the local universe, essentially all isolated 
galaxies with m sta r < 10 9 M® show active star forma¬ 
tion. This is also the case in our models. Therefore, 
at low stellar masses we should be able to measure 
the “native” star-forming main sequence (SFMS), un¬ 
affected by internal quenching processes. We cannot 
reliably reach such low mass limits, but we use galax¬ 
ies with stellar masses between 10 9 and 10 9 ' 5 M© to 
measure the baseline SFMS (we restrict our sample in 
the models to central galaxies for reasons mentioned 
in Section 4.3). We then find the mean log(sSFR) of 
galaxies in this mass range in time bins, tracking the 
evolution of the sSFR of typical star forming galaxies 
across cosmic time. Once this evolution is known, we 
calculate the main sequence slope by measuring the 
change in the mean log(sSFR) between stellar masses 
of 10 9 and 10 10 M©. In a given redshift bin, we use 
the mean low-mass sSFR and derived slope to define 
a mass-dependent main sequence line. We then define 
quiescent galaxies as having less than 25% of the sSFR 
of the main sequence line. Our quiescence divisor for 
a given redshift and stellar mass is given by 


sSFR(z, M„) = 0.25[10 MS(z)+6(i ° 9(M * /M ® )-9 ' 25) ] (2) 

where b is the slope we derived and MS is 
the mean log(sSFR) measured in the low mass bin 
(10 9 M q <M*<10 9 ' 5 M q ). The values of these quanti¬ 
ties are determined separately for the model galax¬ 
ies and for the observed galaxies. The coefficients for 
MS(z) and the values of b in each of our redshift bins 
are listed in Tables 1 and 2, respectively. 

Figure [l] shows the division of star forming and 
quiescent galaxies for the model (including disk in¬ 
stability) and from the observations in the redshift 
bins z ~ 0.06, 0.5 < z < 1.0 and 2.2 < z < 2.6. 
Although the galaxy sample we will use for the re¬ 
mainder of this work includes only galaxies with 
log(M*/M Q )>10.0, here we plot galaxies down to stel¬ 
lar masses of 10 9 M©, since these are the galaxies from 
which our dividing lines are derived. The green line in- 
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dicates our split between star forming and quiescent 
galaxies. The red line is drawn at sSFR = 1/(3£h), 
where tu is the Hubble time at the median redshift of 
the bin, for comparison with alternate dividing lines 
commonly used in the literature. We apply the cut 
derived for the DI model to the noDI model as well, 
since we would like to see how the disk instability 
affects the sSFRs of galaxies within the model and 
that information would be lost if we allowed the cut 
to move between the two models. It is worth noting, 
however, that MS(z) and b are very similar between 
the two models. We can see in all three bins that the 
dividing line has a different slope and normalization 
for the models than for the observations. 


3.1.2 Evolution of the Quiescent Fraction 

Figure [2] shows the evolution of the quiescent fraction 
of galaxies with redshift for galaxies from the SAM, 
both the DI and noDI models, and from observational 
data taken from GAMA and CANDELS. We compute 
1 -cr uncertainties due to field-to-field variance and un¬ 
certainty in observed galaxy properties (stellar mass, 
Sersic index and star formation rate) as follows. Our 
lightcones are about nine times larger than the CAN- 
DELS fields that we are comparing with, so we select a 
model sample from a subsection of the lightcone that 
has comparable area. If we select different CANDELS- 
sized areas from our lightcone to do our analysis, we 
get a measure of the effect of cosmic variance. We also 
calculate the 1-cr error in the quiescent fraction due to 
uncertainties in the estimates of galaxy properties in 
the observational sample. We use quoted uncertainties 
in Sersic index, assume an uncertainty of 0.25 dex for 
star formation rates and use the redshift-dependent 


stellar mass uncertainty of Behroozi et al. (2013b) 


The separate uncertainties due to cosmic variance and 
parameter estimation can be seen in the top panel. We 
add the uncertainty due to each in quadrature and 
apply them to the observations. In the lowest redshift 
bin, the error estimates reflect only the uncertainties 
due to errors in the physical parameters; these uncer¬ 
tainties dominate over the cosmic variance due to the 
large volume probed by GAMA. We note here that 
we are still likely underestimating uncertainties due 
to systematics such as the assumed star formation 
histories of CANDELS galaxies, possible variations in 
the IMF, etc. 

The quiescent fraction of galaxies in the model 
is relatively insensitive (changing by <10% in all red¬ 
shift bins) to the inclusion of disk instabilities in our 
models; as we will see, the net effect of the disk insta¬ 
bility is mainly to create more bulge-dominated galax¬ 
ies. This is due in part to the fact that our disk insta¬ 
bility prescription does not affect gas and limits the 
amount of low-level AGN feedback that is triggered 
by disk instabilities. Both models agree well with ob¬ 
servations at low redshift; for z < 1.2, the fractions 
differ by no more than 0.05 — 0.1. Above this red¬ 
shift, however, the fractions begin to differ by about 
0 .2, with the model predicting fewer quiescent galax¬ 
ies than are observed. Overall, the model exhibits a 
steeper evolution than the observed galaxies, predict- 
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Figure 1. The distribution of observed and model galaxies (log(M*/M0)>9.O) in the plane of stellar mass and specific 
star formation rate (sSFR) in the redshift bins z ~ 0.06, 0.5 < 2 ; < 1.0 and 2.2 < 2 ; < 2.6. The greyscale shows the 
population density with contours overplotted in black. The green line shows our adopted dividing line between star forming 
and quiescent galaxies. In practice, the dividing line is calculated for each galaxy individually based on its stellar mass and 
redshift; the green line is a least mean squares fit to the stellar masses and threshold sSFRs of each galaxy. The red line is 
the 1/3£h dividing line sometimes used in the literature. It is clear that the normalization and slope for the model SFMS 
is different from those for the observations, necessitating the use of a different dividing line. Left panel: Galaxies from the 
SAM. Right panel: Galaxies from the GAMA or CANDELS survey. 
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Table 1. Coeffients for MS(z) = a, 3 t 3 (z) + CL 2 t 2 (z) + ait(z) + ao, the mean log(sSFR) of the main sequence, where t(z) is 
the age of the universe in Gyrs at the redshift of interest. 


Dataset 

a 3 

a 2 

a>i 

clo 

SAM 

- 0.0012 

0.039 

-0.499 

-7.640 

GAMA & CANDELS 

-0.0017 

0.039 

-0.398 

-7.513 


Table 2. Slope derived as described in the text for each of our redshift bins. 


Redshift SAM b GAMA & CANDELS b 


0.006 

< 

z 

< 0.12 

0.021 

-0.303 

0.5 

< 

z 

< 

1.0 

-0.105 

-0.400 

1.0 

< 

z 

< 

1.4 

-0.054 

-0.144 

1.4 

< 

z 

< 

1.8 

-0.241 

-0.130 

1.8 

< 

z 

< 

2.2 

-0.377 

-0.236 

2.2 

< 

z 

< 

2.6 

-0.408 

-0.256 

2.6 

< 

z 

< 

3.0 

-0.487 

-0.370 



Redshift [z] 

Figure 2. The evolution of the quiescent fraction of galax¬ 
ies (log(M*/M 0 )>lO.O) with redshift. The top panel is 
the predicted l-cr uncertainty due to sample variance (red, 
dashed) and due to uncertainty in galaxy parameter esti¬ 
mation (black, solid). These are added together in quadra¬ 
ture and shown plotted on the observational measurements 
in the bottom panel. In the bottom panel the red dashed 
line with squares corresponds to the model including disk 
instabilities, the blue dotted line with triangles to the 
model without disk instabilities, and the black solid line 
with circles to the observations. This convention is used 
throughout this work. Field-to-field variance is expected to 
be negligible in the lowest redshift bin, so here the plotted 
error is entirely due to uncertainties in galaxy properties. 
Overall, the agreement between the model predictions and 
observational results is quite good. Below z ~ 1 . 2 , the frac¬ 
tions predicted by the model differ from the observational 
fraction by no more than 0 . 1 , although at z > 1.2 they 
begin to differ by as much as 0 . 2 , with the model predict¬ 
ing almost no quiescent galaxies. The predicted quiescent 
fraction is affected very little by the inclusion of disk in¬ 
stabilities in the model. 


ing basically no quiescent galaxies at z ~ 3. It seems 
that the model is not quenching galaxies early enough. 
We will discuss possible reasons for this discrepancy 
later. 


In Figure [ 3 ] we examine the mass dependence 
of the quiescent fraction evolution. The behavior is 
similar in each mass bin to the overall behavior in 
Figure [2] At redshifts above ~ 1.2 the model pre¬ 
dicts a smaller quiescent fraction than is observed for 
all three mass bins. This discrepancy gets worse as 
the stellar mass increases. In the highest mass bin, 
the fractions can differ by as much as 40 %. This 
is an extension of the overall high redshift discrep¬ 
ancy in Figure [2] since our model predicts no quies¬ 
cent galaxies (in any mass range) and in general more 
massive galaxies are likely to be considered quiescent, 
the gulf between our model quiescent fraction and the 
observed fraction widens as the masses considered be¬ 
come larger. In the two lower mass bins, the DI model 
actually overproduces quiescent galaxies by as much 
as 0.15 at z < 1.2, but it is in better agreement (within 
0.05) with observations in the highest mass bin. As 
expected, the quiescent fraction increases for galax¬ 
ies with higher stellar mass for both the model and 
the observations. The model also captures the steeper 
evolution of the quiescent fraction for higher masses, 
although again, the evolution in the model is steeper 
than observations for all mass bins. 


3.2 Spheroid-Dominated Fraction 


We now split galaxies into spheroid-dominated and 
disk-dominated populations. We define spheroid- 
dominated galaxies as having Sersic indices greater 
than 2.5, the average of a pure disk (n = 1) and pure 
bulge (n — 4), as has been done in many other studies 


( Shen et al.|2003 Lange et al.|2015 Bruce et al.|2014| 


Mortlock et al.| 2015). We discuss later how making 
this spheroid-domination cut less stringent (at n = 2) 
affects our results. We have also done the same analy¬ 
sis by dividing galaxies at a bulge-to-total mass ratio 
of 0.5. These results are very similar and can be found 
in Appendix B. Figure [4] shows the evolution with red¬ 
shift of the spheroid-dominated fraction of galaxies. 
Here, we see the main effect of the disk instability. 
The noDI model severely underpredicts the fraction 
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Figure 3. Each panel is like Figure[2] but for different bins in stellar mass. Line types and colors are the same as Figure |2] 
The quiescent fraction increases with stellar mass for both the models and observations, with disk instabilities contributing 
more of the quenched galaxies at higher masses. The discrepancy between the model predictions and observations is larger 
in the two higher stellar mass bins. 



Figure 4. The evolution of the spheroid-dominated frac¬ 
tion of galaxies with redshift. Error bars are the 1 — a un¬ 
certainties due to sample variance and uncertainties in ob¬ 
served galaxy properties added in quadrature, as in Fig. [2] 
The separate contributions are plotted in the top panel. 
The model in which spheroids form only via mergers un¬ 
derproduces the fraction of spheroid dominated galaxies at 
z < 2 and does not reproduce the build-up of the spheroid- 
dominated population seen in the observations. The model 
with additional spheroid growth via disk instabilities (DI 
model) is qualitatively in fairly good agreement with obser¬ 
vations (the two agreeing to within ~ 0.1 at all redshifts), 
though the predicted evolution in spheroid fraction is still 
a bit too shallow, with the model overpredicting spheroid- 
dominated galaxies at high redshift and underpredicting 
them at low redshift. 


of spheroid-dominated galaxies at almost all redshifts, 
with the disagreement becoming worse towards lower 
redshifts. The DI model does a much better job of 
matching the observed spheroid-dominated fraction, 
increasing our prediction by almost a factor of two at 
low redshift; within the uncertainties, the fractions do 
not disagree by more than 0.1 at any redshift. How¬ 
ever, the evolution of the spheroid-dominated fraction 
in the model is somewhat shallower than in observa¬ 


tions, so we overpredict the spheroid-dominated frac¬ 
tion at z > 1 and underpredict it at z < 1. 

We once again investigate the mass dependence 
in Figure [5] Again, as expected, at larger masses, 
the spheroid-dominated fraction is greater at almost 
all redshifts, and the spheroid-dominated fraction in¬ 
creases more rapidly with redshift for massive galax¬ 
ies. There appears to be a population of massive, disk- 
dominated galaxies at high redshift in both the obser¬ 
vations and the model. The behavior in all three mass 
bins is reminiscent of the overall behavior in Figure 
[4] except for the high mass, high redshift case. As in 
Figure [4] the disk instability brings the model mainly 
into better agreement, although in the two lower mass 
bins at z > 1 the fractions can differ by as much as 
0.15 — 0.2, which is slightly more than for the overall 
population. 


3.3 Comparison with Previous Results 


We take a moment here to compare with previous 
work that has probed the evolution of the quiescent 
and spheroid-dominated fractions of galaxies. Bram-| 


mer et al. (2011) and Muzzin et al. (2013b) both 


examine the dependence of the quiescent fraction 
on stellar mass across a range of redshifts. |Bram-| 
mer et al. (2011) examine galaxies from the NEW- 


FIRM Medium-Band Survey (NMBS) (|van Dokkum 
et al. |2009l Whitaker et al.|2010 ) while [Muzzin et aL 


( |2013b| ) observe a sample of galaxies in the COS- 
MOS/UltraVISTA field ( Muzzin et al.|[2013a ) over a 
similar redshift range. Our observational results are in 
good agreement with both studies. The quiescent frac¬ 
tion is higher at larger stellar masses and lower red¬ 
shifts and the quiescent fraction of high stellar mass 
galaxies increases more steeply with redshift than that 
of low mass galaxies. The quiescent fraction evolution 
of our high mass bin (10 n M© < M* < 10 11-5 M©) 
is in good agreement with |Brammer et al. ( 2011| ; in 
both cases, the fraction increases from ~ 0.5 at z ~ 2 
to ~ 0.8-0.9 at z ~ 0.5. Muzzin et al. (2013b) is in 
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Figure 5. Same as Figure [3] but for the spheroid-dominated fraction. 


good agreement with our high mass quiescent frac¬ 
tion as well, but also investigates the quiescent frac¬ 
tion down to lower stellar mass so we can compare 
our lower mass bins. These also agree very well. For 
galaxies with 10 lo M© < M* < 10 lo ' 5 M©, the quies¬ 
cent fraction increases from ~ 0 . 2 at 2 :~ 2 to~ 0.4 
at z ~ 0.2-0.5. Meanwhile, the quiescent fraction of 
galaxies with 10 lo ' 5 M© < M* < 10 11 Mq increases 
from ~ 0.4 at z ~ 2 to ~ 0.5-0 .6 at z ~ 0.2-0.5. 

In terms of spheroid-dominated fraction, 


Buitrago et al. (2013) examines the fraction of 
spheroid-dominated galaxies with masses >1O 11 M 0 
from z ~ 3 to the present day. They cover this 
range by combining several different surveys: SDSS 
DR7 flAbazajian et "ah] [20091, POWIR /DEEP2 
(Bundy et al. 2006 Conselice et al.||2007| and the 


GOODS NICMOS Survey (GNS: Conselice et al. 


2011). They find a steady increase in the fraction 


of spheroid-dominated galaxies (n > 2.5) whereas, 
when the same mass cut is applied (as can be seen 
in the right panel of Figure [5]), we predict a sharper 
increase in spheroid-dominated fraction from ~ 20 % 
at z ~ 2.5 to rsj 60% at z ~ 1.5. Because of this, 
we predict a somewhat larger spheroid-dominated 
fractio n than theirs between z ~ 1 and 2 . IBrucel 
et al. (2014) observe galaxies with M*>lO n M 0 in 


the COSMOS and UDS fields from the CANDELS 


survey over the redshift range l<z<3. They use 
bulge-disk decompositions to sort galaxies by B/T 
and compute the spheroid-dominated fraction. They 
find a spheroid-dominated fraction of ~ 0.6 for 

z rsj 1.5 and a fraction of ~ 0.45 for z ~ 2.5. This is in 
very good agreement with our results (again cutting 
at 10 11 Mo and now defining spheroid-dominated 
as having B/T> 0.5); we find a spheroid-dominated 
fraction of ~ 0.45 for z ~ 2.5 and ~ 0.65 for z ~ 1.5. 


3.4 Dividing into Quadrants 

Having divided the sSFR-Sersic plane in halves, we 
now further divide the plane into four quadrants to 
examine the evolution of the populations in each one: 
SFDs, QSs, SFSs and QDs. Figure [ 6 ] shows an exam¬ 
ple of the division of galaxies into quadrants for both 


the model and the observations in the redshift bins 
z ~ 0.06, 0.5 < z < 1.0 and 2.2 < z < 2.6. The 
star formation division line is a least mean squares 
fit to the individual star formation thresholds for 
each galaxy in the redshift bin according to its stel¬ 
lar mass and specific redshift. Figures [ 7 ] and [ 8 ] show 
the distributions of Sersic index and sSFR for our DI 
models and for the observations in the redshift bin 
0.5 < z < 1.0. In both cases, the distributions are sim¬ 
ilar, but not exactly the same. Our model has trouble 
reproducing the strong observed bimodality in both 
quantities; our disk instability creates many galaxies 
of intermediate Sersic index. As we move toward lower 
redshift, the differences between the distributions of 
model and observed galaxies become more significant. 
We will return to this point in the discussion. 

Figure [9] shows the evolution of the fraction of 
all galaxies in each quadrant with redshift for the DI 
model, the noDI model and the observations. We see 
here again the reason for the difference in how the qui¬ 
escent and spheroid-dominated fractions change with 
the disk instability: the disk instability decreases the 
fraction of QDs while increasing the fraction of QSs, 
leaving the quiescent fraction relatively unchanged. 
The two spheroid-dominated populations, however, 
are both increased, leading to the large change in the 
overall spheroid-dominated fraction. The DI model re¬ 
produces the evolution of SFDs to within a few per¬ 
cent as their numbers dwindle due to various transfor¬ 
mative processes. The noDI model predicts too many 
SFDs, the model and observed fractions differing by 
as much as 0.2. The DI model reproduces the observed 
fraction of QSs at z ~ 0.1, but slightly underproduces 
QSs at higher redshifts, although the fractions do not 
differ by more than ~ 0.1. The noDI model under¬ 
produces QSs at all redshifts, to an even larger de¬ 
gree. Both models underpredict the fraction of QDs 
at z > 1.5 (which is again an extension of the over¬ 
all issue seen in Figure [5]) and overpredict them at 
z < 1.5. Once again, the disagreement is worse when 
the model with no disk instability is considered. While 
the models match the observed fraction of SFSs to a 
few percent at redshifts < 1 , they predict too many 
at high redshift, in some cases by a factor of two. It 
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Sersic Index [n] 


Figure 6. The distribution of galaxies in the sSFR-n plane in the redshift bins z ~ 0.06, 0.5 < z < 1.0 and 2.2 < z < 2.6. 
Left panel: Galaxies from the SAM. Right panel: Galaxies from GAMA and CANDELS. The greyscale shows the population 
density in the sSFR-n plane, with contours in black overplotted. The green lines are the dividing lines used in this work to 
identify the four “quadrants” (see text). 


becomes clear when comparing the two models that 
the disk instability is mostly responsible for the ex¬ 
cess of SFSs that we predict at redshifts z ~ 1.5 — 2.5. 
At redshifts higher than this, mergers seem to be¬ 
come increasingly important as a channel for bulge 
growth. We expect some SFSs in the universe to have 


disturbed morphologies due to the process responsi¬ 
ble for making them an SFS. It is possible that some 
of the CANDELS galaxies that would be classified as 
SFSs are dust obscured and are either not detected or 
not considered star forming, leading to an underesti¬ 
mate of the fraction of SFSs. While we do include the 
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Figure 7. Distribution of Sersic indices for model galaxies and CANDELS galaxies in the redshift bin 0.5<z<1.0. Left 
panel: All galaxies. Middle Panel: Star forming galaxies. Right panel: Quiescent galaxies. 



Figure 8. Same as previous figure, but now showing the distribution of sSFR. Left panel: All galaxies. Middle Panel: 
Disk-dominated galaxies. Right panel: Spheroid-dominated galaxies. 


effect of dust extinction in our model, as well as make 
the same iL-band magnitude cut as is used for CAN- 
DELS, the possibility remains that we are underesti¬ 
mating dust extinction. This would cause objects that 
are missed in CANDELS due to the iL-band magni¬ 
tude limit to be included in our model catalogs. It is 
not unreasonable that we would be underestimating 
the effects of dust in these objects in particular, as 
our prescription is based on an undisturbed disk ge¬ 
ometry and does not account for the possibly heavily- 
obscured starbursting systems we are concerned with 
in the SFS quadrant. 

We also note here that changing our cut in Ser¬ 
sic index from n = 2.5 to n = 2, which still distin¬ 
guishes systems with significant bulge components, 
does change our results somewhat as the distribution 
of n in the models is different from the observed dis¬ 
tribution (as seen in Figures [6] and [7]) . The spheroid- 
dominated fraction increases more for the DI model 
than for the observations, especially at higher red- 
shifts. The noDI model is changed very little. When 
looking at different mass bins as in Figure [5] the 
change of the spheroid-dominated fraction of the DI 
model relative to the observations is more pronounced 
in the two lower mass bins than in the highest one. 
Figure [To] shows the evolution of the fraction of all 
galaxies in each quadrant for the morphology cut at 
n — 2. The DI model now underpredicts the fraction 
of SFDs by about 0.1-0.15 at all redshifts and overpre¬ 


dicts SFSs at high redshift by an even larger amount 
(as much as 0.4). The fraction of QSs matches the 
observational results well at z < 1.5 but still under¬ 
predicts these objects at higher redshifts. However, 
qualitatively the results are very similar, so we con¬ 
tinue to use our n — 2.5 cut for the rest of this work. 

Now, knowing both where our model succeeds 
and fails in matching the buildup of these populations, 
we can dig into the model to see which mechanisms 
are responsible for moving our simulated galaxies in 
the sSFR-n plane. 


4 DISCUSSION 

The SAM can provide us with details about galaxy 
formation histories which we cannot glean directly 
from observations; we now examine the statistics of 
events that drive galaxy tranformation and quench¬ 
ing (mergers and disk instabilities) in our models, 
and provide representative examples of how individ¬ 
ual galaxies trace out their histories in the sSFR-n 
plane. 

Figure [Tl] shows density contours for galaxies 
from the SAM (0.5 < z < 1.0) in the sSFR-n plane. 
Overlaid arrows show how different physical processes 
might move galaxies in this diagram. SFDs may merge 
with each other or suffer disk instabilities to form 
bulge-dominated galaxies which then undergo gas de- 
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Figure 9. The evolution of the fraction of galaxies in each quadrant of the sSFR-Sersic plane with redshift. Top left: Star 
forming disk-dominated galaxies (SFD). Top right: Star forming spheroid-dominated galaxies (SFS). Bottom left: Quies¬ 
cent disk-dominated galaxies (QD). Bottom right: Quiescent spheroid-dominated galaxies (QS). Our models qualitatively 
reproduce the trends of a decreasing fraction of SFD galaxies and the increasing fraction of QS galaxies with cosmic time, 
with the DI model in general producing better agreement with the observations. Our models do less well at reproducing 
the observed trends for SFS and QD, predicting mild decrease and increase in these populations, respectively, with cosmic 
time, while in the observations their fractions are nearly constant from 3 > z > 0.1. 


pletion by AGN feedback, leading to quenching of 
star formation. In the SAM this occurs over relatively 
short time scales of several hundred million years. 
Meanwhile, other SFDs may passively evolve, deplet¬ 
ing their gas reservoirs over much longer time scales 
of a Gyr or more, eventually becoming QDs, which 
may then experience dry (gas-poor) mergers which 
puff them up and form QSs. Quiescent galaxies may 
then accrete new gas which allows regrowth of a disk 
component. We now examine the importance of some 
of these processes in the SAM in a bit more detail. 


4.1 How Recently Have Different Types Been 
Disturbed ? 

We would like to know how galaxies in different quad¬ 
rants in the sSFR-Sersic plane are formed or evolve 
to their current state and with the SAM we can di¬ 
rectly measure the time since traumatic events such 
as mergers and disk instabilities. In the top left panel 
we look at the fraction of >10 10 M© 


12 


of Figure 

galaxies of each type which have undergone a recent 
merger, with “recently” being defined as within three 
dynamical times (where tdyn = 27rrdi S k/^disk)- We see 
that SFSs are more likely to have experienced a re¬ 
cent merger at all redshifts, while very few QDs have 


© 0000 RAS, MNRAS 000 , 000-000 












Quenching and Morphological Transformation 17 



Figure 10. Same as Figure[9] but now with a morphology cut at n = 2. 


undergone a recent merger. SFDs and QSs fall in be¬ 
tween. All types of galaxies are more likely to have 
experienced a recent merger at higher redshift. In the 
top right panel, we see that almost no galaxies at any 
redshift have avoided ever having a merger in their 
lifetime. 

In the middle row, we restrict our attention to 
major mergers with a mass ratio >1:3 (the mass 
used to calculate this ratio is the combined bary- 
onic (cold gas+stellar) and dark matter mass within 
2 halo scale radii; see S08) and see that basically no 
QDs have undergone recent major mergers, while the 
fraction of SFDs with recent major mergers is only 
slightly higher (no more than ~ 5%). This is because 
it is very unlikely to experience a major merger and 
still retain enough of a disk to be considered disk- 
dominated within a dynamical time of the merger. 
QSs are slightly more likely to have undergone a re¬ 
cent major merger than SFDs, while 65-80% (depend¬ 


ing on redshift) have undergone a non-recent major 
merger. This is simply because very soon after a major 
merger, the merger-triggered starburst would cause 
the spheroid-dominated galaxy to be classified as an 
SFS. After some time has passed and star formation 
has been quenched, it would be classified as a QS. 
For this same reason, SFSs are the least likely to have 
undergone a non-recent major merger; if the major 
merger wasn’t recent, they’re unlikely to still be star 
forming. SFSs are still most likely to have had a recent 
major merger and that likelihood increases somewhat 
towards higher redshift. The disk-dominated classes 
are actually more likely to have had a non-recent 
major merger than never to have had one at all, af¬ 
ter which they must have regrown a substantial disk 
component. Unsurprisingly, the disk-dominated pop¬ 
ulations are more likely to have never had a major 
merger than the spheroid-dominated populations. 

In the bottom row, we see that the fraction of all 
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Figure 11. Schematic representation of how different physical processes might cause galaxies to migrate in the sSFR-n 
plane. The density distribution of galaxies in the SAM (DI model) for 0.5 < z < 1.0 is shown with contours. 


types which have undergone a recent disk instability 
peaks at z ~ 1.5. We note that this is in qualitative 
agreement with the peak of the clumpy fraction of 
galaxies (in the mass range 9.8<log(M*/M©)<10.6) 
found in Guo et al. (2015). We cannot make a quan¬ 


titative comparison with these results as we have no 
way to estimate clumpiness in our models, but clumpy 
galaxies are expected to be associated with minor 
mergers and disk instabilities (D ekel et al.|| 2009). 


Our star forming classes are more likely than our 
quiescent classes to have undergone a recent disk in¬ 
stability, but the fraction of SFSs plummets towards 
higher redshift, signaling the fact that mergers seem 
to be the dominant bulge-growing channel at high red- 
shift, with disk instabilities increasing in importance 
as mergers become less frequent, as has been spec¬ 


ulated in previous works (Parry et al. 2009 De Lu 


cia et ah||2011 Fontanot et al.||2012 ). QSs are much 

more likely to have never had a disk instability, or to 
have had one longer ago (since it has likely been a 
while since they had a disk). QDs are just as likely to 
have had a disk instability recently, not recently, or 
not at all, suggesting that disk instabilities don’t play 
a huge role in their evolution. Finally, a significant 
fraction of all types (40-60%) have never experienced 
a disk instability. This fraction increases steeply to¬ 
wards high redshift, presumably because the merger 
rate increases and it is less likely for a galaxy to go 


undisturbed long enough to develop an instability on 
its own. 

We note that there is a large amount of uncer¬ 
tainty in the highest redshift bin for our quiescent 
classes. There are very few galaxies classified as qui¬ 
escent and so we suffer from small number statistics 
in that bin. 


4.2 Individual Galaxy Histories 


To illustrate how individual galaxies evolve, we now 
inspect the evolutionary tracks of four galaxies se¬ 
lected from the SAM, which end up in the four differ¬ 
ent quadrants of the sSFR-n plane at z — 0. Here 
we use bulge-to-total mass ratio as our proxy for 
morphology, since the tracks in the sSFR-morphology 
plane are much easier to see this way (and because the 
results are very similar: see Appendix B). Figure 13 
shows the evolutionary path of a galaxy with a fairly 
quiet history that ends up as an SFD. This galaxy 
has a mass of ~ 10 lo ' 8 M© at z — 0. The top panel is 
its track in the sSFR-morphology plane, color coded 
by the age of the universe. The bottom and middle 
panels are the evolution with time of the morphology 
and sSFR respectively. We see that this galaxy has a 
few mergers early on after which its evolution is en¬ 
tirely due to accretion of new material, allowing it to 
continue forming stars. The decrease in B/T follow- 
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Figure 12. Top row: Left panel: Fraction of model galaxies in each quadrant which have undergone a recent (<3£dyn) 
merger with mass ratio >1:10. Middle panel: Fraction which have undergone a merger on a timescale >3tdyn- Right panel: 
Fraction which have never undergone a merger. Middle row: Same as top row, but for major mergers (>1:3). Bottom Row: 
Same as top and middle rows, but for disk instabilities. The fraction of galaxies that have suffered a recent merger declines 
with cosmic time from z ^ 2 to the present, while the fraction that have experienced a recent disk instability peaks at 
around z ~ 1.5-2. 


ing the merger events is due to the regrowth of a disk 
component. 

Figure [14] is a somewhat striking example of a 
galaxy being pummeled repeatedly by mergers until 
it is almost entirely bulge dominated, after which it 
finds itself unable to form more stars because of the 
black hole it has grown over the course of its traumatic 
history; the black hole is now keeping any remaining 
gas too hot for star formation through radio-mode 
feedback. This can also be effected by one (or more) 
big major merger(s). In both cases, it is likely that the 
system will make an appearance as an SFS for a time 
before quickly evolving into the QS quadrant. Once a 
galaxy falls into this quadrant, it tends to stay there, 


except for very rarely when it collides with a gas-rich 
galaxy, at which point it might briefly return to the 
SFS quadrant before quickly using up all of its new 
gas and falling back down again. This can be seen in 
the very short spikes of star formation accompany¬ 
ing the last two mergers in the bottom panel. This 
galaxy is far more massive than the SFD above, with 
a stellar mass of ~ 1O 11 ' 8 M 0 , which is due to all of 
the merger events it has experienced and is consistent 
with many of the most massive galaxies in the uni¬ 
verse being QSs. We note that the QS population in 
the models is much more concentrated toward lower 
Sersic index than that from the observations, as can 
be seen in Figure [6] Our model is still not producing 
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Redshift [z] 



Figure 13. Evolution of sSFR and B/T mass ratio for a 
galaxy which is classified as a star forming disk-dominated 
galaxy with a mass of ~ 1O 1CL8 M0 at z = 0. Top panel: 
Evolutionary track in the sSFR-B/T plane, color coded by 
age of the universe. Minor mergers (<1:3), major merg¬ 
ers (>1:3) and DIs are indicated by diamonds, circles and 
triangles, respectively. Middle panel: Evolution of B/T 
mass ratio with time. The red dash-dotted line indicates 
a major merger and the blue dotted lines indicate minor 
mergers. The solid red line is our division between disk- 
dominated and spheroid-dominated. Bottom panel: Evo¬ 
lution of sSFR with time. The solid red line is our division 
between star forming and quiescent. The galaxy remains 
disk dominated, due to its quiet accretion history, and the 
SFR gradually declines due to the declining cosmological 
accretion rate. 


enough very bulge-dominated galaxies. As the most 
massive bulges are believed to be the result of merg¬ 
ers, this may indicate that we are still underestimating 
the role of merging. 

Figure [15] shows the evolution of a galaxy that 
ends up as an SFS. Its final mass is ~ 10 11 M©, which 
is more massive than the SFD considered above be¬ 
cause of its more active merger history. We see here 
what a major merger can do in terms of bulge growth: 
the first major merger this galaxy experiences gives 
it a substantial bulge component. After each bulge 
growth episode, the galaxy begins to regrow a disk, 
causing B/T to decrease steadily. This is the case 
for many galaxies in the SAM, as long as they have 
the gas to form new stars or continue to accrete new 
gas from the IGM. Changes in sSFR between merger 
events for this system are largely due to the interplay 
between “normal” star formation and new gas accre¬ 
tion. If this galaxy had not undergone a major merger 
recently, it would not have been considered an SFS, 
as the steady regrowth of its disk would have caused 
it to be classified as disk-dominated instead. 

The SFS quadrant is more of a way station than 
a destination in galaxy evolution, with galaxies either 
quickly evolving back towards the SFD quadrant by 
regrowing their disks, or evolving downwards into the 
QS quadrant as they quench (this system is very close 
to being classified as a QS). This can be seen in Fig¬ 
ure [6] where the SFS quadrant appears more to be 
made up of the tails of populations in the SFD and 
QS quadrants than to be a distinct population. This 
is something upon which the models and observations 
appear to agree, and is an explanation for why the 
evolution of the overall spheroid-dominated fraction is 
less steep than for the quiescent fraction: there is high 
turnover in the SFS quadrant so the fraction in that 
state is relatively constant, meaning that the buildup 
of the spheroid-dominated population relies mainly on 
the steady buildup of QSs. Meanwhile the quiescent 
fraction is built up by the steady growth of both the 
QS and QD populations. This interpretation seems 
to be corroborated by the CANDELS-based study of 
Rizer et al. (in prep.), in which visual morphologi¬ 
cal classifications are used. They find that while their 
QS population builds up steadily, their SFS popula¬ 
tion remains relatively constant, suggesting that bulge 
growth is leading to star formation quenching in many 
cases. 

Figure [16] represents one possible path to becom¬ 
ing a quiescent disk-dominated galaxy. Some quies¬ 
cent disks in our model are galaxies with quiet histo¬ 
ries such as the one seen in Figure [T3| above, but which 
have run out of gas (perhaps because they are more 
massive and have a harder time accreting new mate¬ 
rial) or can’t form stars with the gas they do have (see 
below for a brief discussion). However, many of our 
disk-dominated quiescent galaxies really aren’t very 
disk-dominated but in fact are systems which have 
had a somewhat more eventful history similar to that 
of a QS. These end up as QDs by one of two ways. In 
the first case, the events which lead to star formation 
quenching and bulge formation do not form enough 
of a bulge for the system to be considered spheroid- 


© 0000 RAS, MNRAS 000 , 000-000 

























Quenching and Morphological Transformation 21 


dominated. However, the galaxy still falls into the 
QD quadrant in the same way we expect SFSs to mi¬ 
grate sometimes to the QS quadrant after a traumatic 
event. In the second case, the galaxy does become an 
SFS, but retains gas to regrow a disk sufficient to be 
considered disk-dominated before falling into quies¬ 
cence. 


We see a mixture of these two fates in Figure 
[Tg] while its bulge growth episode was not enough to 
make it spheroid-dominated, the event which caused 
the bulge growth was enough to cut this galaxy off 
from new gas, causing it to become quiescent. We 
also see substantial disk growth in its decreasing B/T. 
This combination of both scenarios causes this galaxy 
to be somewhat rare: a QD which is both very disk- 
dominated and very quiescent. This galaxy is very 
massive at z — 0 with a stellar mass of ~ 10 11-7 M©. 
QDs are often quite massive due either to their merger 
histories or prolonged star formation which has led 
the galaxy to be unable to accrete new gas due to 
virial shock heating. In general, the average B/T for 
QDs is only slightly higher than for SFDs. However, 
as QDs with lower star formation rates are considered, 
the average B/T becomes larger, as these systems are 
the result of the QS-like paths described above. The 
very disk-dominated QDs are likely to have relatively 
higher sSFRs. This does bring up a larger point about 
dividing lines in general: some galaxies which fall be¬ 
low our dividing line are still forming stars, albeit at a 
slower rate than the majority of galaxies of their mass. 
They are “quenched” in the sense that their star for¬ 
mation rate is lower than expected, but they are not 
truly “quiescent” as is the case for some of our galax¬ 
ies, which have had their star formation completely 
turned off. There is also the possibility, mentioned 
above, that some of our quiescent galaxies might be¬ 
gin to form stars again due to a wet merger or new gas 
accretion. In the case of a wet merger, this is likely to 
be short-lived, but in the case of new gas accretion, 
it can lead to a whole new life for a galaxy. Having 
clarified that, however, we believe this is happening 
in both the models and the observations, so it should 
not bias our results. 


It appears that the QDs are much like the SFSs 
in that the population is a combination of tails of the 
SFD and QS populations. Curiously, the gas fractions 
of QDs in the SAMs are not systematically smaller 
than those of SFDs of the same mass. It appears that 
the quenching leading to very disk-dominated QDs 
has two possible origins, as mentioned above: a low 
gas accretion rate, or an extended low surface-density 
gas disk which is inefficient at making stars. In the 
models considered here, only gas that is above a crit¬ 
ical surface density is allowed to participate in star 
formation; galaxies with larger than average angular 
momentum form more extended disks, which have a 
larger fraction of their gas sitting below this critical 
surface density. There is observational evidence in the 
local universe for these gas rich but relatively quies¬ 
cent disks ( Lemonias et al.||2014 Schiminovich et al, 


2010). Thus while the SFSs may be understood as a 


transient population, a step in the path from SFD to 
QS, it seems that QDs may be a static population. 


In addition, our model QDs on average have 
smaller stellar and black hole masses than our QSs. 
They are also more likely to be satellites than QSs. 
In our lowest redshift bin, about 75% of model QSs 
are central galaxies, while only 50% of QDs are. The 
very disk-dominated QDs ( B/T < 0.1) don’t exist in 
our model before z ~ 2. When compared with the 
rest of the QDs, they have even smaller black hole 
masses, as they likely have not undergone any events 
which would have triggered AGN feeding. They are 
also even more likely to be satellites; in our lowest 
redshift bin, only 25% of very disk-dominated QDs 
are central galaxies. We note that the likelihood of be¬ 
ing a satellite galaxy is not the only reason that QDs 
tend to have smaller stellar and black hole masses; 
this trend is observed even when only central galaxies 
are considered. 

What these analyses and evolutionary tracks 
show is that the transformative processes which affect 
galaxies take them all over the map (and the sSFR-n 
plane). It is likely too simplistic to tell a simple story 
about two star forming disk galaxies colliding or one 
of them buckling under its own weight and triggering 
feedback which produces a nice, dead elliptical galaxy. 
These processes (mergers, disk instabilities, accretion 
of new gas) likely work together, sometimes in tandem 
and sometimes at cross purposes. It appears that a 
complex history with multiple transformative events 
is the norm rather than the exception and galaxy his¬ 
tories don’t necessarily look like the arrows in Figure 

m 


4.3 How Might a Different Implementation 
of Disk Instability Change Our Results ? 

As mentioned above, our treatment of disk instabil¬ 
ities is based on rather dated simulations of isolated 
disk galaxies which are not in a cosmological context, 
and may not capture all of the relevant physics. Some 
of the questions associated with disk instabilities in¬ 
clude: 1) What is the most relevant criterion for deter¬ 
mining the onset of a disk instability? 2) What hap¬ 
pens to the gas and stars in the disk when it becomes 
unstable? 3) How efficiently do disk instabilities feed 
a nuclear black hole? 

As an aside, we have so far elided over a pos¬ 
sibly important distinction. There are two kinds of 
physical mechanisms that are commonly referred to 
as “disk instabilities” in the literature, although in 
one case this is something of a misnomer, as we ex¬ 
plain below. In “violent disk instabilities” (VDI; Dekel 
et al. 2009), the disk becomes globally unstable, lead 
ing to the formation of clumps of stars and gas, which 
may migrate to the center of the galaxy, building the 
spheroid. As the giant clumps orbit within the disk, 
even if they disrupt before reaching the center, they 
may drive inflows of gas into the galaxy nucleus, via 
the same sort of physics as merger-induced nuclear in¬ 
flows, again leading to growth of the spheroid through 
in situ star formation, and feeding of the black hole 
( Bournaud et al.||2011 ). The second kind of “disk in¬ 
stability” involves the secular transfer of angular mo¬ 
mentum outwards, and mass inwards, again leading to 
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the building of a central compact and dynamically hot 
structure, and is accompanied by the formation of a 
bar (Kormendy & Kennicutt 2004). The term “insta¬ 
bility” is a misnomer here, as the disk essentially re¬ 
mains in dynamical equilibrium. The relatively crude 
morphological statistics used in this analysis are not 
able to distinguish between edge-on bars and bulges, 
so some of the ‘spheroids’ we count in the observa¬ 
tions may actually be bars. In addition, the nuclear 
structures that are formed via this secular process do 
not have the same properties as “classical” bulges, and 
are sometimes called “pseudobulges”. See Kormendy 
&; Kennicutt (2004) for a detailed discussion of the 


differences between classical bulges and pseudobulges; 
the most germane for our purposes here is that pseu¬ 
dobulges do not obey the same scaling relationships 
with the SMBH mass as classical bulges (Kormendy 
fcHo| 2013]), suggesting that black hole feeding and/or 
feedback may operate differently. 

One major limitation of our approach is that we 
based the disk instability criterion on the properties 
of the stellar disk only, and only stars are moved from 
the disk to the spheroid when the disk is deemed un¬ 
stable. This is not what is seen in modern cosmological 
hydrodynamic simulations, in which VDI are ubiqui¬ 
tous at high redshift ( Ceverino et al.|2010| Mandelker 


et al. 2014), and as mentioned above are associated 


with strong nuclear inflows of gas as well as stars. In 
addition, VDI can lead to significant quenching even 
in the absence of associated AGN feedback (Gammie 


2001 Dekel V Burkert 2014 Forbes et al. 2014), which 


is not accounted for in our current models. Perhaps 
this could lead to a higher fraction of quiescent galax¬ 
ies at high redshift, in accord with observations. Addi¬ 
tional complications are the possible stablizing effects 
of a pre-existing bulge or central mass concentration, 
and possible triggering of disk instabilities by minor 
mergers. 

P14 attempted to model a scenario closer to the 
modern VDI picture with their “Stars+Gas DI” 
model. They found that the results were very simi¬ 
lar to the “Stars DI” model which is why we con¬ 
sider only that model here (this implementation of 
DI is similar to the one that is most commonly ap¬ 
plied in other SAMs in the literature). In addition, the 
P14 “Stars+Gas DI” model is still very arbitrary 
and simplified. In future work, we plan to carry out 
a thorough analysis of state-of-the-art hydrodynamic 
simulations to develop a more detailed and physical 
treatment of disk instabilities in SAMs. 


4.4 Other Possible Model Improvements 

Aside from our treatment of disk instabilities, there 
are several ways in which we could improve our model 
in order to better capture what we believe is occuring 
in real galaxies. Our implementation of satellite strip¬ 
ping does not account for any change in morphology 
of satellites and many which are not destroyed are 
stripped of their gas and remain disk dominated. Only 
central galaxies are supplied with new gas in our ac¬ 
cretion and cooling model, so this may artificially in¬ 
crease our fraction of quiescent disk-dominated galax¬ 


ies. While this effect dominates at low stellar masses 
and should not be as important in the mass range 
we consider in this work, M* > 10 lo M©, it is worth 
noting (and is the reason we neglected satellites when 
determining the typical model main sequence star for¬ 
mation rate above). We do point out that the same 
analysis done for the model excluding satellite galax¬ 
ies leads to very similar results; the largest difference 
is in the low redshift QD fraction, which in the lowest 
bin decreases from ~ 25% to ~ 17 — 18%. 

While our merger prescription is based on nu¬ 
merical hydrodynamic simulations of binary mergers 
(Hopkins et al. 2009b), we treat all mergers as dis¬ 
crete events, while in the real universe there may be a 
complex interplay between merger events if the galaxy 


has not had time to relax between them (Moster et al. 


2014). This could have an effect on both the star for¬ 


mation rates and the morphologies of our post-merger 
remnants, especially of those that have had somewhat 
active merger histories. 

In addition, in our current models galaxies are 
primarily kept quiescent by the “maintenance mode” 
type radio mode feedback, which becomes important 
only at relatively late times. Winds driven by radia- 
tively efficient accretion (“bright mode” feedback) are 
assumed to be able to remove cold gas from the ISM 
but have no effect on the hot gas surrounding galaxies. 
However, recent cosmological simulations have shown 
that momentum-driven winds associated with bright 
mode accretion actually modify the hot gas profile and 


significantly retard cooling over long timescales (Choi 


et al. 2014 2015). Including this physics might also 


help us to produce more quenched galaxies at high 
redshift, and could also suppress the re-formation of 
disks via cooling. 

Finally, there is the issue of putting the models 
and observations on the same footing when it comes 
to quantitative comparison of morphologies; we would 
like to improve the way we assign morphology to our 
model galaxies, namely by using the masses (or lumi¬ 
nosities) and sizes of our disk and bulge components 
to generate mock images which can then be processed 
like real observations (including noise, point-spread 
function, etc) and assigned morphological classifica¬ 
tions. This way each galaxy would have its profile 
measured separately and other effects such as incli¬ 
nation angle could be taken into account. Comparing 
the observed and predicted fractions of morphologi¬ 
cally disturbed galaxies would also be interesting, but 
obviously requires us to be able to quantify this in the 
SAM in a way that can be compared with observa¬ 
tions. This may be possible using our merger statis¬ 
tics in combination with a library of numerical simu¬ 
lations, as in Lotz et al. (2011); see also recent work 
by Snyder et al. (2014). 


5 SUMMARY AND CONCLUSIONS 

We have studied the coevolution of star formation 
rate and morphology from z ~ 3 to the present by 
examining the buildup of galaxies in the four quad¬ 
rants of the sSFR versus Sersic index (n) plane. We 
have compared galaxies with stellar mass >10 lo M© 
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from the “Santa Cruz” semi-analytic model outlined 
in S08, S12 and P14 with galaxies observed as part of 
the GAMA and CANDELS surveys. Our conclusions 
are as follows: 

• Our models qualitatively reproduce the increas¬ 
ing fraction of quiescent galaxies since z ~ 2 seen 
in observations, and produce excellent quantitative 
agreement with observations at z < 1.2. At higher 
redshift, the models underproduce the fraction of qui¬ 
escent galaxies relative to observations. 

• Our model in which spheroids are built solely 
through mergers (noDI) predicts an evolution in the 
fraction of spheroid dominated galaxies that is much 
too mild compared with observations. This model also 
underproduces spheroid-dominated galaxies at z < 2 
compared with observations. Adding a channel for 
bulge growth via disk instabilities (DI model) leads 
to much better agreement with the observed evolu¬ 
tion of the spheroid-dominated fraction, although still 
produces slightly flatter evolution than observed. 

• The quiescent fraction is largely unaffected by the 
inclusion of disk instabilities, with changes <10%. We 
note, however, that our current disk instability model 
may be underestimating the change in quiescent frac¬ 
tion, especially at high redshift. 

• Our models further qualitatively reproduce the 
observed evolutionary behavior of four classes of 
galaxies defined by both star formation activity 
and morphology: star forming disk-dominated, qui¬ 
escent spheroid-dominated, star forming spheroid- 
dominated and quiescent disk-dominated. In both the 
observations and in our models, the fraction of star 
forming disks decreases over time while the fraction 
of quiescent spheroids increases. In the observations, 
the fractions of both star forming spheroids and qui¬ 
escent disks remain nearly constant from z ~ 3-0. 
Models predict a stronger decrease in star forming 
spheroids and a stronger increase in quiescent disks 
with redshift than is seen in the observations, but the 
predicted fractions are not off from the observed frac¬ 
tions by more than ~ 0.2 and in most cases are off by 
less. 

• In our models, star forming-disk dominated 
galaxies are galaxies which have had very quiet his¬ 
tories. They have avoided major mergers and if they 
have experienced any merger or disk instability ac¬ 
tivity, they have recovered by accreting new gas and 
regrowing a disk. 

• Again in our models, quiescent spheroid- 
dominated galaxies are likely to have either undergone 
at least one extreme major merger or many smaller 
mergers or disk instabilities. In either case they have 
built up a substantial bulge component and AGN 
feedback has made it impossible for them to accrete 
significant amounts of new gas, eventually leading to 
cessation of star formation. 

• Star forming spheroid-dominated galaxies seem 
to be a short-lived population. Truly spheroid- 
dominated star forming galaxies are indicators of a 
recent trauma, as they are still experiencing a post¬ 
trauma starburst. At this point, they can regrow a 
disk with their remaining gas reservoir or through the 


accretion of new gas. In the absence of new gas, they 
can deplete their gas reservoirs and become quiescent. 

• Quiescent disk-dominated galaxies are a combi¬ 
nation of two populations: disk-dominated galaxies 
which have stopped accreting gas (in some cases due 
to environmental effects) and galaxies with extended 
low-surface density gas disks, which are inefficient at 
forming stars. 

Despite the room for possible improvements to 
our model described above, the qualitative similarity 
between the buildup of our model populations with 
those of observed galaxies gives us confidence that we 
are beginning to capture the complicated interplay of 
several processes which lead to the diversity of galax¬ 
ies and their evolution over time. 
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APPENDIX A: CONVERSION FROM B/T 
TO SERSIC INDEX 


As part of this work, we compare the morphologies 
of model galaxies to observed galaxies. To do this, 
we convert bulge-to-total stellar mass ratios (B/T), 
which our model naturally outputs, to Sersic indices 
from single component fits, which are directly compa¬ 
rable to observations. To do this, we use a lookup table 
generated from synthetic galaxies which are made up 
of an exponential (n — 1) disk and a bulge with n— 4. 
A Sersic index and effective radius are derived by fit¬ 
ting to the two-component profile for a wide range 
of B/T and rbuige/^disk (which we hereafter refer to 
as “rbd”). The lookup table takes in B/T and rbd 
and outputs a Sersic index and effective radius for 
the composite system. Since the values are discrete, 
we have interpolated between the table values to gen¬ 
erate our Sersic indices. Appendix A of Lang et al.| 
(2014) illustrates that the relationship between B/T 
and n derived from these synthetic galaxies matches 
well with the relationship derived from CANDELS 
galaxies with 2-component bulge+disk fits. Here we 
present the mapping between B/T and Sersic index 
in order to illustrate the relationship between the two. 
Then, in the next section, we will show that the results 
of the analysis presented in the main text are largely 
unchanged when done in terms of B/T rather than 
Sersic index. We also refer the reader to Figure A6 
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in Lang et al. (2014), where this analysis is done for 


the observed galaxies with bulge-disk decompositions. 
They carry out their analysis in terms of both bulge- 
to-total stellar mass ratio and H-band light ratio. We 
have used the lookup table generated in terms of the 
mass ratio, but we would expect results in terms of 
light ratio to be qualitatively similar (see Appendix B 
of Lang et al.|2014 ). 

To test the SAM predictions for the bulge and 
disk sizes, we use the r-band bulge+disk decompo¬ 
sitions of SDSS galaxies performed by |Simard et al.| 
(2011). We trimmed the original catalog of 1,123,718 


galaxies down to 618,186 galaxies by applying the fol¬ 
lowing selection cuts: 0.005 < z < 0.12, 0.0 < B/T < 
1.0, 0.5 < Tlpure < 8.0, Mr > —99, M r , e rr > —99, 
Mr, pure ^ 99, Mr, pure, err > 99, r*bulge,eff ^ 0, 
t disk,eff > 0, and r pU re,eff > 0, where the subscript 
“pure" refers to single-component (pure) Sersic fits 
(which were also computed for the galaxies). The 
bulge+disk decompositions were fit simultaneously 
in the r- and g-band in order to minimize errors, 
and the assumed model was a de Vaucouleurs bulge 
(^buige = 4) with a pure exponential disk. The fits 
were done using the Galaxy Image 2D (GIM2D) pro¬ 
gram; see |Simard et al. (2002) and Sim ard et al.| ( |20lT ) 
for further details about the fitting procedure and 
outputs. In general, the model predictions and ob¬ 
servational results are similar, except in the lowest 
B/T bin, where our model predicts more compact 
bulges relative to the disk sizes than is seen in the 
observations. This excess is seen to a lesser degree in 
the other B/T bins as well. However, as the discrep¬ 
ancy is the largest for disk-dominated galaxies (where 
the radial size of the bulge component will have little 
impact on our results), we conclude that our model 
should produce reasonable predictions for the compos¬ 
ite Sersic indices in most cases. In the future, it will 
be interesting to compare the SAM results with the 
sizes and B/T ratios obtained from multi-component 
bulge-disk decompositions, which are starting to be¬ 
come available. 

Figure [TS] shows the best fit curves to the re¬ 
lationship between B/T and n in bins of rbd. If 
rbd< 0.4, for a given B/T a larger rbd will lead to a 
lower Sersic index (unless the galaxy also has very low 
B/T, in which case n is mostly concentrated between 
~ 0.5 — 1.0 anyway). Above an rbd of ~ 0.4, there is 
a nearly one-to-one mapping between B/T and Sersic 
index. However, as we saw in Figure [TT| in the SAMs 
and in nearby galaxies, most galaxies have rbd< 1. 
We see here that a split in B/T instead of n will lead 
to the selection of slightly different sets of galaxies 
because the bulge radius to disk radius ratio causes a 
spread in Sersic index for a given bulge-to-total mass 
ratio and vice versa. 


tions are carried out, and may be more easily com¬ 
pared with predictions from other theoretical mod¬ 
els. Our dividing line between star forming and quies¬ 
cent remains the same, but our condition for spheroid 
domination is now B/T > 0.5. Figure 


19 


shows the 

evolution of the spheroid-dominated fraction of galax¬ 
ies as in Figure [4] The observed spheroid-dominated 
fraction is still derived using a Sersic index n — 2.5 
to make the cut. We can see that the evolution of 
model galaxies is similar to the evolution in Figure [4] 
However, splitting by B/T we predict more spheroid- 
dominated galaxies for z < 2 and fewer for z ~ 3. 
The predicted fraction does not vary by more than 
0.1 at any redshift. This variation suggests that there 
are more galaxies at high redshift with smaller values 
of rbd than at low redshift. 

The evolution of the quadrant fractions in Fig¬ 
ure [20] is extremely similar to that in Figure [5] The 
fraction predicted in each of the spheroid-dominated 
quadrants is slightly larger than in the Sersic index 
case. Figure [2l] is like Figure [12] The only significant 
difference when splitting by B/T instead of n is in the 
disk instability plot. More SFSs and QSs have had re¬ 
cent disk instabilities as defined by B/T than by n. 
This is because the change in definition is not likely 
to affect galaxies that are very clearly disk-dominated 
or spheroid-dominated. The galaxies that have inter¬ 
mediate Sersic indices or B/T ~ 0.5 are the ones that 
are traded back and forth depending on definition (be¬ 
cause of their rbd values), and a large fraction of these 
are created by the disk instability. This is also why 
the results for the DI model seem to be affected more 
strongly by morphological definition than the noDI 
model. 


APPENDIX B: RESULTS USING B/T 

In this section we present the main (morphology- 
dependent) results of the above analysis again, this 
time using B/T as our morphological parameter. 
These results may be able to be compared with fu¬ 
ture observational analyses, if bulge-disk decomposi- 
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Figure 14. Same as previous figure, but for a galaxy which 
is classified as a quiescent spheroid-dominated galaxy at 
z = 0 with a mass of ~ 10 1 1-8 Mq. Although this galaxy 
has suffered one major merger at very high redshift (z = 
5.5), the build-up of its dominant spheroid occurs through 
a sequence of multiple minor mergers. This is quite typi¬ 
cal. The build-up of the spheroid is accompanied by growth 
of the SMBH, leading to strong radio-mode feedback that 
shuts down cooling and, eventually, quenching of star for¬ 
mation. 


Figure 15. Same as previous figure, but for a galaxy which 
is classified as a star forming spheroid-dominated galaxy 
at z = 0 with a stellar mass of ~ 10 11 Mq. This galaxy 
has had an extremely active history, with multiple major 
mergers, several minor mergers, and a disk instability. A 
recent major merger triggered a strong burst of star forma¬ 
tion. This galaxy would probably appear morphologically 
disturbed. 
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Figure 16. Same as previous figure, but for a galaxy which 
is classified as a quiescent disk-dominated galaxy at z = 0 
with a stellar mass of ~ 1 O 11 , 7 M 0 . 
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Figure 17. Distribution of rbd values for SAM galaxies with 0.06 < z < 0.12 and galaxies from SDSS with bulge-disk 
decompositions from Simard et al. | ([2011) in bins of B/T. 
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Figure 18. Fit to the relationship between B/T and Sersic index in bins of bulge radius/disk radius. For compact bulges, 
the Sersic index is a function of both B/T and rbd. As rbd increases, the relationship becomes degenerate. 
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Figure 19. The evolution of the spheroid-dominated frac¬ 
tion of galaxies with redshift, now with spheroid domi¬ 
nation defined as B/T > 0.5. The observed galaxies are 
still split by Sersic index at n = 2.5. Error bars are the 
1 — a uncertainties due to sample variance in the mod¬ 
els and uncertainty in observed galaxy properties added in 
quadrature. The separate contributions are plotted in the 
top panel. 
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Figure 20. Same as Figure [ 9 ] but now with model galaxies split by B/T = 0.5. Observed galaxies are still split by Sersic 
index at n = 2.5. Top left: Star forming disk-dominated galaxies. Top right: Star forming spheroid-dominated galaxies. 
Bottom left: Quiescent disk-dominated galaxies. Bottom right: Quiescent spheroid-dominated galaxies. 
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Redshift [z] 

Figure 21. Same as Figure [12] but with a morphology cut at B/T = 0.5. Top row: Left panel: Fraction of galaxies in 
each quadrant which have undergone a recent (<3£dyn) merger with mass ratio >1:10. Middle panel: Fraction which have 
undergone a merger on a timescale >3£dyn- Right panel: Fraction which have never undergone a merger. Fractions are now 
determined using the B/T cut. Middle row: Same as top row, but for major mergers with mass ratio >1:3. Bottom row: 
Same as top and middle rows but for disk instabilities. 
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